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Some Results of a Multiple Ship Survey of the Gulf Stream” 


By PCR GLISmARe andl. eve WORTHINGTON, 
Woods Hole Oceanographic Institution 


(Manuscript received 18 December, 1950) 


Introduction 


The first multiple ship survey of the Gulf 
Stream Area between Cape Hatteras and the 
Grand Banks of Newfoundland took place 
during the period from the 6th to the 23rd o 
June 1950. The survey was given the code 
name “Operation Cabot” by the U. S. Navy 
Hydrographic Office which acted as head- 
quarters during the planning and operational 
phases of the survey. The civilian scientist in 
charge of the technical phases of the program 
was Dr RicHarp H. FLEMING, Director, 
Division of Oceanography, U. S. Navy 
Hydrographic Office. 

The purpose of Operation Cabot were: 


(a) To obtain a synoptic plot of the Gulf 
Stream’s path between Cape Hatteras and 
the Grand Banks; 

(b) to determine the rate of change in the 
position of the Gulf Stream; 

(c) to measure the surface velocities of the 
Gulf Stream; 

(d) to observe what effect the Gulf Stream 
has on the conditions in the lower 
atmosphere; 

(e) to observe the formation and some of the 
history of the eddies presumably formed 
when a larger “meander” or “wave” 
in the Gulf Stream is cut off; 

(f) to obtain observations of 


Gulf 
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Stream from a plane and from ships 
simultaneously. 


The following groups took part in the 
survey. 


U. S. Navy Hydrographic Office, Washing- 
tons DAC 
U. S.S. SAN PABLO (AGS-30); Captain, 
A. J. PETERSON Comdr. USN.; Scientist 
in Charge, Dr R. H. FLEMING. This vessel 
acted as flagship for the fleet and Dr FLEMING 
was task unit chief scientist. 
VUSSESZREEOBOME (AGS-so); Cap- 
tain, R. V.R, Basser, Ir, Comdr. USN; 
Scientist in Charge, ROBERT ABEL. Two 
PB4Y-2 Aircraft; Pilots Ensign F. M. 
GLAZIER, USNR and Lt J. G. STEWART, 
USNR;; Scientific Observer, W. V. Kier- 
HORN. 
Naval Research Establishment, Halifax, Nova 
Scotia. 
H.M.C.S. NEW LISKEARD; Captain, 
„Et W., Wo-Maceoll. REN. Scentise ın 
Charge, Dr Wırrıam L. Forp. 
U. S. Fish and Wildlife Service, Washington, 
DEG. 
ALBATROSS III; Captain, JoHN COLLINS; 
Scientist in Charge, V. B. CoLTon. 
Woods Hole Oceanographic Institution, 
Woods Hole, Mass. 
R/V ATLANTIS; Captain, ADRIAN LANE; 
Scientist in Charge, MARTIN POLLAK. 
R/V CARYN; Captain, JoHN Pixs; Scientist 
in Charge, Dr B. H. Kercuum. 
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Atlantic Oceanographic Group, St. Andrews, 

New Brunswick. 
The original plans called for the participa- 
tion of the CNAV WHITETHROAT. 
Because of last minute difficulties the 
WHITETHROAT was unable to join in 
the survey. Dr H. B. HACHEY and Dr 
L. M. LAUzIER of the organization helped 
to plan the operation and Dr LAUZIER was 
on the scientific staff aboard the NEW 
LISKEARD. 

Scripps Institution of Oceanography, La Jolla, 
Calif. 

Dr Watter Munk of this organization 
helped plan the operation and was on the 
Scientific Staff of the flagship. 

U. S. Navy Office of Naval Research, Wash- 

ington, D.C. 
Lt Comdr. R. L. Danrorr and J. K. Knauss 
of this organization helped plan the opera- 
tion. Lt Comdr. Dautorr was chief meteoro- 
logist ot the task unit. 

Although not an integral part of Operation 
Cabot the work done by the U.S. Coast 
Guard Research Vessel EVERGREEN, while 
on routine Ice Patrol for the International 
Ice Patrol, contributed valuable data. Dr 
FLoyp SOULE was Scientist in Charge >n 
the EVERGREEN. 


Operation Cabot was a splendid example of 
cooperative effort and credit for the success 
of the survey must be distributed among all 
the many scientists and seamen who took part. 


Operational plans 


The planners of Operation Cabot were 
faced with a new type of oceanographic 
investigation. Not only were six ships and a 
plane involved but new instruments for ob- 
taining data with ships underway were to 
be used. A very large amount of data was to 
be collected and, although much of it would 
be collected automatically by seamen tech- 
nicians, the scientific staffs would be kept 
busy supervising and evaluating the results. 

Oceanographic surveys have traditionally 
been made according to a pre-arranged plan, 
and in an area such as the one under considera- 
tion, the measurements have been planned 
along “sections” or at points on a “grid”. 
Recent surveys of the Gulf Stream (IsELIN and 
FUGLISTER, 1948) have indicated that such a 
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method is inherently wasteful of time and 
effort where answers to specific questions 
are sought, because many stations do not 
provide pertinent information. Nevertheless, 
detailed plans and instructions were formulated 
for Operation Cabot. This was done primarily 
because it was felt that without such pre- 
arranged plans the success of the survey 
would be entirely dependent on successful 
communications. Considering the variety of 
ships and equipment involved dependable 
communications were by no means certain. 
Also the planners realized that if communica- 
tions were successful the flagship could alter 
the plans during the operation as might be 
desirable. 

Briefly the plans for the survey were as 
follows. The first phase was to be a general 
reconnoitering of the Gulf Stream from Cape 


. Hatteras to the Grand Banks. The ships would 


be deployed approximately 150 miles apart 
along the Gulf Stream, zig-zagging into and 
out of the current and underway continuously. 
Bathythermograph, Loran, and weather data 
were to be obtained at half-hour intervals. 
Surface - current measurements were to be 
made every hour with the Geomagnetic 
Electrokinetograph (von ARX 1950), hereafter 
referred to as the GEK. Positions of the Stream 
were to be transmitted every six hours to 
the flagship. 

The second phase was to consist of simultane- 
ous “sections of hydrographic stations along 
prearranged meridians of longitude. Bathy- 
thermograph, Loran, and weather observations 
were also to be made along the “sections”. 

Phase three was to be a detailed study of a 
cestricted area determined by the results of 
the first two phases. 

Throughout the operation continuous sonic 
depth records were to be made, and a track 
chart maintained with the depth records re- 
ferenced to the track chart by hourly markings. 

Fortunately communications were very 
nearly perfect during the entire operation. 
The first phase of the operation went off as 
scheduled and the scientific staff aboard the 
flagship was able to plot the position of the 
Gulf Stream from off Cape Hatteras to 50° west 
longitude. For the remainder of the Opera- 
tion the flagship directed the movements of 
all vessels and the pre-arranged plans were in a 


large part discarded. 
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An enormous amount of data were collected 
during the 17 days of the survey and the present 
paper will cover only a portion of the results. 


Definition of terms 


Early in the planning and operational stages 
of Operation Cabot it became evident that 
precise definitions were needed for the various 
terms used in association with the Gulf Stream. 
The frequent references to the “cold wall”, 
“edge of the Stream”, “warm core” and 
“front” led to a certain amount of confusion 
and misunderstanding. The term “inner edge” 
was most frequently used and most variously 
interpreted. This confusion is caused primarily 
because, although the words “Gulf Stream” 
denote a current, they also imply a distinct 
water mass and secondarily because water 
masses that may be motionless are included 
as part of the Stream because they lie below 
a surface current. 

Since the Gulf Stream is a boundary or 
front in the western North Atlantic between 
the slope water and the Sargasso Sea we 
may define it as follows; it is a continuous 
band stretching from the continental shelf off 
Cape Hatteras to the soth meridian of longi- 
tude, south of the Grand Banks of Newfound- 
land. This band consists of a pronounced 
pressure gradient between the warm highly 
saline water to the south and the colder 
fresher water to the north. Using this defini- 
tion then, the inner and outer limits or edges 
of the Gulf Stream can be defined as the 
points where this pressure gradient becomes 
zero. These points can be located only if 
deep, closely spaced temperature and salinity 
data are obtained and the cross current pressure 
gradients calculated. Also, because of the 
large eddies found both north and south of 
the Stream, any section made across the area 
must be long enough to ascertain whether or 
not more than one pronounced pressure 
gradient exists. If only one is found it defines 
the Gulf Stream, but if more than one are 
located then the position of the Stream cannot 
be determined by that single section. 

Not to be confused with the inner or left 
hand edge of the Stream is the temperature- 
salinity boundary at the surface. This generally 
abrupt change that occurs to the left of the 
“warm core” may or may not coincide with 


the left hand edge of the Gulf Stream as 
defined above. This applies also to the color 
boundary and the long thick lines of Sargas- 
sum frequently seen on the surface; all of these 
surface phenomena are apparently associated 
with shear zones to the left of the “warm 
core” but they are not necessarily coincident 
with the left hand or inner edge of the Gulf 
Stream. 

The “warm core” is defined here as that 
part of the Gulf Stream containing water 
warmer than the water at the same depth to 
the right, facing down stream, of the current. 
This “warm core” is generally 300 to 400 
meters deep with the maximum temperature 
anomalies at a depth of about 100 meters. 

The word “front” is considered synonymous 
with the pronounced pressure gradient and 
therefore with the Gulf Stream itself. 

The term “cold wall” dates back to 1845 
and is still frequently used to denote the “inner 
edge of the Stream” or, according to CHURCH 
(1937), “the temperature gradient between 
the slope water and the Gulf Stream”. Accord- 
ing to these definitions it could equally well 
be called the “warm wall” though in neither 
case do we have anything resembling a wall. 
This temperature gradient exists at different 
depths across the entire width of the Gulf 
Stream and therefore cannot be considered 
as something separate or adjoining the Stream. 
Because of the misleading connotations of the 
term “cold wall” it will not be used in this 


paper. 


Synoptic charts 


On this survey no attempt was made to 
locate, precisely, the inner and outer edges 
of the Gulf Stream. The area covered is shown 
in fig. 1. It is obvious that to obtain closely 
spaced, deep hydrographic stations over this 
large area, in the limited time available, was 
impossible even with six ships. However a 
study of recent surveys show that the limits 
of the Stream can be determined with a fair 
degree of accuracy from the temperature 
structure in the upper 200 meter layer. A 
temperature-depth profile across the Stream 
shows the following characteristics; facing 
down stream the left hand edge of the Stream 
coincides with the point where the temperature 
at 200 meters reaches a minimum or levels off. 


F. C. FUGLISTER AND 


L. V. WORTHINGTON _ 


75° 


60° 


x -13000 FOOT BATHY THERMOGRAPH 
e -HYDROGRAPHIC STATIONS 


N 


| | 
x1 ORIN, 


o 1 Past PA. 
a0- i nr 7 
xl ah, OX. HYDROGRAPHIC 
(IX) ET D DRE STATIONS 
Api SRE ESS OCCUPIED BY THE 
4 


USCGC EVERGREEN 
ON ICE PATROL 
9-20 JUNE, 1950 


so —_—_— 1 


70° 65° 


55° 50° wnoi SEPT 1950 


Fig. 1. Allships’ tracks on Operation Cabot, 6th to 23rd of June 1950. 


The abrupt temperature boundary at the 
surface may or may not be directly over this 
point. To the right the isotherms slope down- 
ward until the “axis of the warm core” is 
reached. This is the point where the mean tem- 
perature of the upper 200 meter layer is a 
maximum. Actually the “warm core” reaches 
to depths of 300 to 400 meters, but since the 
maximum anomalies are found at about 100 
meters, it has proven convenient to use the 
mean of the upper 200 meter layer to show 
the position of this core. To the right of the 
“warm core” the isotherms in the surface 
layer rise until the temperature at a depth of 
200 meters is again a minimum or levels off. 
This point is taken as the right hand or outer 
edge of the Stream. The horizontal tempera- 
ture and salinity gradients at all levels are 
quite gradual on this side of the Stream so 
that this “edge” is not nearly so well defined 
as the “inner edge”. 

In order to obtain synoptic charts of the 
Gulf Stream the first step was to calculate, 
for each bathythermogram, the mean tempera- 
ture of the upper 200 meter layer. 3,521 
bathythermograph observations were obtained 


during this survey and temperatures were read 
at eight points on 3,000 of these. Since the 
bathythermograph grids were in the Fahren- 
heit scale the labor involved in converting to 
the Centigrade scale was prohibitive. This 
explains the various scales used in this paper. 

The positions and mean temperature (de- 
grees Fahrenheit) of all bathythermograph 
observations were plotted geographically on 
daily charts and the results contoured. The 
mean temperature isotherms for the 8th of 
June, over the western portion of the Gulf 
Stream are shown in fig. The solid 
lines in this figure show the ships’ tracks and 
give an idea of the amount of interpolation 
that was generally necessary in drawing the 
curves. In relatively few of the 165 “sections” 
made was the entire width of the Gulf Stream 
covered, but in all cases the abrupt gradient to 
the left of the “warm core” was crossed. The 
sections that did cross the “warm core” 
indicate that it did not vary markedly in width 
throughout its entire length from Cape 
Hatteras to the soth meridian. Fig. 3 shows 
the positions of the “warm core” during the 
first and last periods of the survey. The most 
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Fig. 2. Mean temperature, °F., in the upper 200 meter layer of the Gulf Stream, on the 8th of June 1950. 


striking feature of this chart is, of course, the 
two large waves and the deep trough separating 
them in the center of the area. After the 12th 
of June the chief emphasis of Operation Cabot 
was to observe the changing conditions in 
and around this trough. As shown in the 
figure, the trough eventually “broke off” 
from the Stream forming, as it did so, an 
elongated cyclonic eddy to the south. 
Although Operation Cabot was the first 
survey to cover the entire distance from off 
Cape Hatteras to so’ west longitude, other 
shorter surveys, show that large waves in the 
Gulf Stream may occur over a wide range of 
longitude and that conditions, such as those 
found on “Cabot” are not exceptional. Fig. 4 
is a composite picture, showing the positions 
of the maximum horizontal temperature gra- 
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Fig. 3. Positions of the “warm core’’ of the Gulf Stream 
during the first and last periods of Operation Cabot. 


dient at a depth of 100 meters, as obtained on 
all Gulf Stream surveys made since 1946. Be- 
sides the numerous large waves, this diagram 
shows another elongated, cyclonic eddy cen- 
tering at 66° west similar to the one found on 
Operation Cabot. 

The variable positions of the major waves 
in the Stream is stressed here because, during 
Operation Cabot, there was a strong feeling 
aboard the flagship that the bottom topo- 
graphy, near 39° north, 61° west, was the 
primary cause of the deep trough that existed 
in this area. Further investigation may show 
that the series of ridges and seamounts in this 
locality always affect the course of the Gulf 
Stream, but the occurrence of large waves 
and eddies in regions where the bottom is 
relatively flat indicates that bottom topo- 
graphy must be of secondary importance in 
the formation of such phenomena. 


Temperatures of the warm core 


Throughout the entire length of the Gulf 
Stream the mean temperature of the upper 
200 meter layer in the warm core was greater 
than 68° F. An unpublished technical report 
by FUGLISTER on the monthly average tempera- 
tures in the region south of the Stream and 
north of Bermuda shows the mean temperature 
for this layer in June to be 67.8° F. Since the 
Stream moves in an easterly direction for a 
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Fig. 4. Positions of the maximum cross-current temperature gradients at a depth of 100 meters, from all surveys, 
1946 to 1950. 


distance of about 1,200 miles it is evident 
that it must transport over this entire distance 
water from more southerly latitudes. 

The mean temperature at the axis of the 
“warm core’ ranged from 08.3 F. to. 935.5 F. 
As can be seen in fig. 2 these temperatures 
did not decrease uniformly down stream but 
what might be called “gobs’’ of warmer water 
showed up at intervals along the Stream. 
These “gobs” were observed as far east as 
the 57th meridian. The more consistent 
temperatures on either side of. the “warm 
Core “indicate “that ‘these cobs Mare not 
caused by internal wave action. If this is true 
then these masses of warmer water, appearing 
at intervals in the Stream, may be evidence of 
a pulsing action in the current which in turn 
may have a bearing on the formation of 
waves and eddies. The pulsing could con- 
ceivably be caused by short period variations 
in the strength of the Trades or in variations 


in the barometric pressure over the Gulf of 
Mexico. If a series’ of temperature sections, 
say two a day for a period of a month, could 
be run across the Straits of Florida some of 
these speculations could be checked. An exami- 
nation of existing bathythermograph sections 
in the Florida current between Key West and 
Cape Canaveral has not revealed any “gobs” 
but it is not possible, in view of the inadequacy 
of the data, to conclude that they do not 
occur there. 


Rate of change in the position 
of the Stream 


In the first period of the survey four ships 
entered the Gulf Stream off Cape Hatteras at 
24 hour intervals. During this interval the 
Stream shifted, at this point, at a rate of 4.5 
miles a day toward the southeast. Numerous 
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Fig. 5. Positions of the Gulf Stream in the western portion 

of the survey at the beginning and end of the operation. 

Mean temperature, °F., of the upper 200 meter layer, 

for the 8th of June solid lines; for the 21st and 22nd of 
June dashed lines. 


sections made near 72° west showed that the 
crest of the first wave in the path of the 
Stream moved toward the east at a rate of 
approximately 11 miles per day. The overall 
change in positions of the western portion of 
the Stream, from the 8th to the 22nd of 
June, is shown in fig. 5. The daily synoptic 
charts, of the mean temperature in the upper 
200 meter layer, show that all the changing 


positions of the Stream can be accounted for. 


by lateral movements of not over 11 miles 
per day. The synoptic charts for the 18th and 
roth of June, for the region around the trough 
(see fig. 6b), may at first glance indicate a 
lateral movement more rapid than 11 miles 
per day if the motion is considered as a north— 
south one. Actually the water on either side 
of the trough converged at a rate less than ıı 
miles per day in cutting off the cyclonic eddy 
to the south. 

Various changes in the Gulf Stream posi- 
tion, are shown in fig. 5, 6a, and 6b. 
Considering the huge volume of the Stream, 
this shifting of its position by as much as 
11 miles per day must produce large vertical 
motions in the water masses adjacent to the 
Stream. Although the evidence for these shifts 
in the Stream is based on observations of the 
upper 200 meter layer, a sufficient number of 
deep stations were made during Cabot and 
other recent Gulf Stream surveys to show 
that the water at least down to 2,000 meters 
is affected. 

Because of the above measurements it is 
now possible to evaluate, with more confi- 
dence, the results of earlier Gulf Stream 
surveys made with single ships. The changes 
in positions of the Stream as found by the 
Canadian Research Vessel NEW LISKEARD 
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Fig. 6 a. Positions of the Gulf Stream, roth to 15th June. 

Shaded areas denote a mean temperature in the upper 

200 meter layer of less than 65° F. Thin line is 68° F. 

isotherm. Small numbers indicate the maximum observed 
values of the warm core. 
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Fig. 6 b. Positions of the Gulf Stream, 16th to 20th 
June. Shaded areas denote a mean temperature in the 
upper 200 meter layer of less than 65° F. Thin line is 
68° F. isotherm. Small numbers indicate the maximum 
observed values of the warm core. 
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Mean temperature, °F., in the upper 200 meter layer on the 17th of June. Current directions from the 


GEK. 


in November 1948 and the large eddy, near 
66° west longitude, surveyed by the ATLAN- 
TIS in June 1947 can be accounted for in 
detail by lateral shifts of the Stream of the 
order of 11 miles per day. It is obvious that 
several ships are required to obtain a synoptic 
picture of any very extensive portion of the 
Gulf Stream but, if it does not shift its position 
at a rate considerably greater than 11 miles 
per day, a single vessel can still be used to 
advantage in studying a limited area. In cither 
case, however, the system of using a pre- 
arranged “grid” of hydrographic stations 
must be avoided. To obtain one section of 
closely spaced, deep stations across the Stream 
requires about 3 days time; a ship moving at 
10 knots could make from 10 to 14 bathy- 
thermograph and GEK sections in this time. 


Surface currents 


Up to this point the positions of the Gulf 
Stream have been deduced from temperature 
observations and it has been assumed, first, 
that the pressure gradients coincide with the 
temperature gradients, and secondly, that the 
currents flow along the isobars. One of the 
most interesting results of Operation Cabot is 
the evidence that the surface current vectors 
and the isotherms (showing the mean tempera- 
ture of the upper 200 meter layer) are parallel 
where the cross current temperature gradient 
is sharpest. 

As mentioned previously, the mean tem- 
perature of the upper 200 meter layer was 
plotted for each day. This was done without 
any reference being made to the observed 
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current directions or speeds. At the completion 
of this work the current vectors, as computed 
by Mr von Arx from the GEK data, were 
superimposed on the isotherms. Examples 
of the results are shown in figures 2, 7, and 8. 
To the left of the warm core where the iso- 
therms are most closely spaced it appeared, 
in nearly all cases, as though the isotherms 
had been drawn along the current vectors. 
Current vectors were obtained for 136 out 
of the 165 sections crossing this temperature 
gradient and the isotherms and vectors were 
parallel in 95 % of the cases. 

To the right of the warm core the direction 
of the current as obtained from the GEK data 
was extremely variable, on many occasions a 
narrow counter current being indicated. In 
one region, within the wave east of the deep 
trough, all sections that crossed the warm core 
showed a counter current. Actually this 
current proved to be part of a cyclonic eddy 
with a cold core similar to, but not as cold as, 
that found in the larger eddy formed by the 
deep trough. It has not been determined 
whether this smaller cyclonic eddy originated 
as a broken off segment of the Gulf Stream 
or was simply a counter current that became 
enclosed or cut off within the large wave to the 
east of the deep trough. 

The current speeds obtained with the GEK 
are being subjected to a close analysis by Mr 
von Arx. He reports that a correction must be 
applied, in the vicinity of the warm core, where 
the depth of the current affects the values 
obtained by his instrument. This correction 
factor has not been determined. The ship’s 
drift data calculated from Loran observations, 
in the western half of the area, showed the 
band of maximum current speed to average 
between 4 and 5 knots. East of 60° the Loran 
data were not sufficiently accurate, to justify 
current calculations by this method. The 
uncorrected GEK readings showed a falling 
off of velocities toward the east, the maximum 
south of the Grand Banks being 85 % of the 
maximum found near Cape Hatteras. In 
general the maximum, uncorrected, velocities 
“were found to the left of the axis of the 
warm core. The maximum shear was also 
observed to the left of the warm core. The 
uncorrected readings gave an average shear 
of 2.5 10 4sec-1 in the steepest part of the 
velocity profile. 
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Fig. 8. Temperature and current pattern of “Edgar” 

during the first and last periods of its observed develop- 

ment. Mean temperature, °F., of the upper 200 meter 

layer. Current directions from the GEK. Speeds in deci- 
meters per second. 


Eddies 


For a number of years it has been noted that 
large cyclonic eddies are frequently encount- 
ered south of the Gulf Stream (IsELIN, 1940, 
ISELIN and FUGLISTER, 1948). In June 1947 
the ATLANTIS circumnavigated one of 
these eddies. Centering at 66° west longitude 
it extended 200 miles east and west and 60 
miles north and south. It was believed that 
these eddies were formed when a large loop 
or trough in the Gulf Stream became cut off 
and separated from the Stream. That this 
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was indeed the case is shown by the observa- 
tions made on Operation. Cabot. 

The day by day changes in the Gulf Stream’s 
course in and around the trough at 61° west 
longitude are shown. in fig. 6a and 6b. 
In. these figures the water with a mean tem- 
perature, in the upper 200 meter layer, of 
less than 65° F. is indicated by shading. The 
heavy line is the 65° F. isotherm, the light 
line the 68° F. isotherm. The light lines are 
the approximate boundaries of the warm 
core. The maximum surface currents were 
observed where these two isotherms (the 
65° and 68°) are shown close together, in 
other words the maximum speeds appeared 
to the left of the warm core. Details of the 
temperature and current pattern are shown 
in fig. 7 and 8. 

For how long a period the trough existed 
before the roth of June is not known but it is 
surprising that, considering the narrowness of 
the “neck”, it took nine days for it to break off. 
During this period the long narrow “neck” 
was meandering on a relatively small scale 
and the temperature characteristics in the 
cold central portion were continuously chang- 
ing. Since there is comparatively little hori- 
zontal flow associated with the cold core these 
temperature changes must indicate vertical 
motion. On three occasions temperature 
profiles across the “neck” indicated that it 
was about to break off but each time subse- 
quent observations showed the recurrence of 
colder water in the upper 200 meter layer. 
The final occlusion in this shallow layer 
occurred on. the roth of June and on the 20th 
a deep hydrographic station showed that the 
break had occurred at all levels. 

Even. before it broke off, the eddy, which 
was christened “Edgar” during the cruise, 
started to turn so that at the end of the survey 
its major axis was oriented east and west. 

The smaller cyclonic eddy northeast of 
“Edgar” was crossed six times. Details of its 
size and shape are lacking but there is enough 
evidence to show that during the period it 
shrunk in size and moved steadily toward the 
northwest. 

The following questions concerning these 
eddies were left unanswered by Operation 
Cabot. Do the eddies once they have broken 
away from the Stream move, and if so in 
what direction? How long do they last? 


Are elongated eddies typical? Do elongated 
eddies eventually break up into smaller, more 
nearly circular masses? Is there any preferred 
point of origin? If it is assumed that the region 
where “Edgar” was formed is such a “pre- 
ferred point of origin”, because of the ridges 
and seamounts, then the detached eddies must 
frequently move ‘in an upstream direction, 
since they have been found as far west as 73° 
west longitude. If on the other hand we assume 
that the detached eddies remain motionless or 
move down stream then it is even more diffi- 
cult to explain how an eddy such as “Edgar” 
could be formed in the relatively confined 
area near Cape Hatteras or south of it. All 
the available data show that the Gulf Stream 
does not start to meander until after it leaves 
the continental shelf at Hatteras. 

An envelope drawn about the curves shown 
in fig. 4 indicates that it is not until the 
Stream reaches at least as far east as the 7oth 
meridian that waves large enough to form an 
eddy could develop. Indeed, such an envelope 
would make it appear that, if waves of the 
size of those found on Operation Cabot are 
necessary to form an eddy, then they do not 
develop until the Stream reaches close to the 
longitude where “Edgar” was formed. 

These considerations combined with what 
is known of the bottom topography indicate 
that the eddies formed by the Gulf Stream 
originate to the east of 63° west longitude and 
after being detached move in an upstream 
direction. 


Possible branching of the Gulf Stream 


Only two ships operated in the region east 
of 57° west longitude during this survey. The 
USCGC EVERGREEN made a bathyther- 
mograph and GEK section along the soth 
meridian from the Grand Banks south to 
38° 30’ north latitude then, returning northward 
she made a series of hydrographic stations as 
well as bathythermograph and GEK observa- 
tions. The HMCS NEW LISKEARD obtained 
bathythermograph and GEK sections into and 
out of the Gulf Stream to 50° 20’ west longi- 
tude. Both of these tracks are clearly seen in 
fisivr. 

The current map as calculated aboard the 
EVERGREEN from the data collected from 
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the 9th to the 2oth of June is shown in fig. 9. 
At about 42° north the EVERGREEN crossed 
a warm current flowing toward the east. This 
current was, in every respect, similar to the 
current that the Ice Patrol finds near this 
latitude each year. This current is referred to 
as the North Atlantic current (Sours, 1938). 
Continuing south the EVERGREEN crossed 
a minor counter current near 40° north and 
then at 39° it entered another easterly stream. 
Three days later the NEW LISKEARD 
showed that the southernmost current was the 
Gulf Stream which she had followed con- 
tinuously from the west. 

In August and again in September 1949 the 
Gulf Stream was located by the ATLANTIS 
at the 39th parallel south of the Grand Banks. 
An unpublished chart of the average depth of 
the 10°C. isotherm in the North Atlantic 
indicates that the Gulf Stream is south of 40° 
at this longitude. On the basis of these data 
and the records of the Ice Patrol it seems 
evident that at 50° west longitude there are 
always two easterly currents of major magni- 
tude, one generally located near 41° north the 
other at about 39° north. According to the 
definition given for the Gulf Stream it is 
quite evident that the southernmost current 
is the Stream. The data from Operation 
Cabot shows clearly that it forms part of a 
continuous current reaching from Cape Hat- 
teras to 50° west longitude. What then is the 
origin of the current to the north and in what 
way are these two current related? Is the one 
to the north a band or a branch of the Gulf 
Stream or is it entirely independent of the 
Stream? 

The temperature-depth profile calculated 
from the bathythermograph and hydrographic 
station data collected by the EVERGREEN is 
shown in fig. 10. The current vectors as 
obtained by the GEK are plotted across the 
top of the profile. The deep observations were 
spaced approximately 30 miles apart while 
the temperatures down to 250 meters were 
obtained at 5 mile intervals. The right hand 
or southern end of the section is a typical 
temperature profile across the Gulf Stream. 
The position of the warm core is evident. The 
narrow weak counter current on the right 
hand side of the warm core is indicated by 
the vectors. The uncorrected current speeds 
(cm/sec) are shown above the vectors, the 


46°) + 


GRAND BANKS 


i 


44° 


40°) i 1 
GENERAL CHART 
COVERING 
ICE PATROL AREA 
ABOUT THE 
GRAND BANKS 
CURRENT MAP 
FROM DATA COLLECTED 9-20 JUNE, 1950 


| 407 


a | 9714 


Re Te 


DYNAMIC TOPOGRAPHY REFERRED 
TO THE 1000 DECIBAR LEVEL 


4175] CIRCLES DENOTE STATION POSITIONS 
| | | 

Er 52° Ey 50° 49° 48° a7" 46° 45° 44° 
WHO! SEPT 1950 


Fig. 9. Current chart obtained from the International 
Ice Patrol. 


maximum speed of 167 cm/sec being to the 
left of the warm core. 

The northern end of the profile shows the 
cold waters of the Labrador current. Between 
the Gulf Stream and the Labrador current 
this profile could be interpreted as showing a 
large anticyclonic eddy that may have been 
formed originally by a wave in the Gulf 
Stream breaking off on its northern. side. 
However the Ice Patrol current maps, Over a 
period of years, have indicated that this 
current near 41° north latitude is of a per- 
manent character and flows generally from 
the northwest and therefore it cannot be part 
of an eddy detached from the Gulf Stream. 

The other possibility is that the northern- 
most current is a band or branch of the stream. 
The authors’ interpretation of the data obtained 
on this survey does not show any such branch- 
ing of the Gulf Stream although, because of 
the relative scarcity of observations east of 
the 58th meridian, it does not preclude the 
possibility. The most simple interpretation is 
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Fig. 10. Temperature, °F., profile from the Grand Banks south on the soth meridian. Based on combined bathy- 
thermograph and hydrographic data. Current directions from GEK and velocities in centimeters per second. See 
fig. 9 for station positions. 


that the northernmost current coincides with 
the pressure gradient between the slope water 
and Labrador current water and the two 
easterly currents at 50° west are independent 
of each other. 


Other observations 


The continuous sonic depth records obtained 
during the survey have been turned over to 
the U.S. Navy Hydrographic Office to be 
incorporated in their bottom topography 
charts. During the cruise the U.S.S. SAN 
PABLO surveyed two seamounts not shown 
on the Hydrographic Office charts. The 
first located at approximately 38° 54’ north 
and 60° 28’ west had a least depth of 725 
fathoms. This was tentatively named San 
Pablo Seamount. The second at approximately 
37° 29’ north and 59° 52’ west was flat-topped 


with an elongated zone some 8 miles long and 
3 miles wide, and less than 700 fathoms. 
The minimum recorded depth was 675 
fathoms. This was tentatively named Cabot 
Seamount. Both seamounts were explored 
by running at least three sections across them 
on different headings. The seamount shown 
on the charts at approximately 39° 00’ north 
and 61° 00’ west was found to have depths 
somewhat less than indicated and that it 
too was flat-topped at a depth of about 500 
fathoms. 

Some excellent colored photographs showing 
shear zones, long thick lines of Sargassum, 
abrupt changes in wave types, clouds and 
color changes associated with the Gulf Stream 
were obtained from the aircraft. The reports 
of the Stream’s positions, based on such 
visual evidence, were generally in close 
agreement with shipboard observations. The 
abrupt changes in the appearance of the sea 
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surface did not however always coincide with 
the inner edge of the Gulf Stream as defined 
in this paper. Although dye marker experi- 
ments carried out by the plane showed clearly 
that the visual boundary was both a velocity 
discontinuity and a line of convergence, 
the ships’ observations showed that there 
may be more than one such surface boundary, 
none of which necessarily coincide with the 
inner edge of the Stream. 

Near 41°00’ north 63° 00’ west, while 
following a long dense line of Sargassum, the 
plane encountered the U.S.S. SAN PABLO 
as it was crossing from south to north through 
the line of weed. The temperature-depth 
profile as obtained by the ship on this section 
is shown in fig. 11. The line of Sargassum 
indicated convergence at the sea surface. 
Visual observations of the ship’s wake and 
GEK data showed that there was an abrupt 
change in velocity at this line but both the 
temperature profile and the GEK showed 
that the ship was still in the current after it 
had passed to the north side of the weed. Fig. 11 
shows that the line of Sargassum was about 
s miles from the abrupt temperature change 
at the surface and about 20 miles from the 
“inner edge” of the Gulf Stream as defined 
by the pressure gradient. 

Color differences associated with the Gulf 
Stream were apparent.to the plane only in the 
westernmost area, west of 71° west longitude. 
East of this point lines of weed, differences in 
wave type, and formations of “typical” Gulf 
Stream cumulus clouds had to be depended 
upon to locate the Stream. Neither the cumulus 
clouds nor the weed were always present and 
the wave patterns, of course, depended on the 
wind velocities and were most pronounced 
with moderate winds blowing against the 
current. In general, visual observations failed 
to identify the Stream east of the 65th me- 
ridian. : 

The authors hope that Mr KIELHORN, the 
observer on the plane, or the Hydrographic 
Office will publish a complete report on the 
"aircraft participation in Operation Cabot 
reproducing some of the excellent photographs 
obtained. 

The large amount of weather data collected 
by the planes and ships during the survey have 
not as yet been studied. These data, the sea 
surface temperature data and more details of 
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Fig. 11. Temperature profile across 
ıoth of June 1950. 


the Gulf Stream, 


the surface current data will form the basis 
of future reports. 


Summary 


Operation Cabot, the first oceanographic 
survey of its kind to date was successful mainly 
because: 

(a) Communications were maintained be- 
tween all ships and the flagship throughout 
the survey, enabling the staff aboard the 
flagship to direct and coordinate the move- 
ments of the various vessels as the synoptic 
picture of the Gulf Stream developed. 

(b) Instruments for taking numerous and 
rapid measurements were available; Loran for 
navigation, the bathythermograph for water 
temperatures, the geomagnetic electrokineto- 
graph for current velocities, and the planes 
for visual observations of the Gulf Stream’s 
positions. 

(c) Highly qualified scientific personnel 
were available to man the six ships and the 
plane. 

In this paper the Gulf Stream has been 
defined as the continuous band, reaching from 
the continental shelf off Cape Hatteras to the 
soth meridian, where a pronounced pressure 
gradient exists between the warm, highly 
saline water to the south and the colder, fresher 
water to the north. The inner and outer edges 
of the Stream are the points where this pressure 
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gradient becomes zero. Abrupt changes in 
temperature and current velocities that occur 
at the surface do not necessarily coincide with 
these boundaries. 

Daily synoptic plots of the Gulf Stream 
show that the current may shift its position at 
a rate of 11 miles per day. Two synoptic 
plots of the Stream are shown, one for the 
first and one for the last periods of the survey. 
The outstanding feature of the first plot is 
the two large waves separated by a deep 
trough near 61° west. The second plot shows 
that the trough became cut off and formed a 
cyclonic eddy south of the Gulf Stream. The 
breaking off of the trough was observed in 
detail and confirmed the theory that the large 
cyclonic eddies frequently found south of 
the Gulf Stream are formed in this manner. 

A smaller, weaker cyclonic eddy was 
observed to the East of the deep trough. 
Although it is assumed that this eddy was 
formed in the same manner as the larger one, 
it is pointed out that it may be a portion of 
the counter current found on the outer side 
of the Stream, that became enclosed and cut 
off within the large wave east of the deep 
trough. 

Eddies are probably always formed when 
the waves in the Gulf Stream reach a certain 
optimum size. It is estimated that the waves 
cannot become large enough to form eddies 
west of longitude 70° and probably not west 
of longitude 63°. There is the further possibility 
that the waves only reach this critical size in 
regions where ridges in the bottom topo- 
graphy affect the course of the Stream. 

There is evidence that the cyclonic eddies 
broken off to the south of the Stream move in 
a westerly or upstream direction. 


The surface currents are shown to flow 
parallel to the isotherms where these isotherms 
are closely spaced in the region to the left of 
the warm core of the Gulf Stream. In the 
right hand portion of this core the current 
directions are variable, frequently showing 
a narrow weak counter current. 

Temperatures of the warm core indicate 
that the Gulf Stream transports water from 
south of Cape Hatteras to the longitude of 
the Grand Banks, the mean temperature of 
the upper 200 meter layer decreasing approxi- 
mately 6° F. over this 1,200 mile distance. 
“Gobs” of warmer water, possibly indicating 
a pulsing action in the Stream, were found at 
intervals in the warm core from Cape Hatteras 
to 57° west longitude. 

To the south of the Grand Banks two 
easterly currents are found and the southern 
one is shown to be the Gulf Stream. No 
evidence was obtained on this survey to show 
that the northerly current was in any way 
connected with the Gulf Stream. 
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Abstract 


In this and a following paper an investigation is made into the factors which are responsible 
for the characteristic “‘streakiness’’ of certain atmospheric and oceanic current systems. As a 
first step a search is made for the factors which might give rise to the striking vertical concentra- 
ton of momentum observed in many well-defined atmospheric and oceanic jet streams. A 
straight, parallel, heteorogeneous and incompressible current with a free surface, hydrostatic 
pressure distribution and a continuous density field is analyzed theoretically to determine the 
particular distribution of velocity with depth which, for a prescribed volume transport within 
each isopycnic layer, leads to a minimum value for the transfer of momentum. This requirement 
leads to a certain “‘critical’’ relationship between the velocity and density distribution which 
in the case of a two-layer system reduces to the well-known critical relative speed at an internal 
fluid boundary. In the general case of a continuous density distribution it is found that currents 
which are subjected to momentum losses and hence tend to approach the “critical” velocity 
distribution in many cases must shrink vertically and develop a sharp velocity maximum at or 
below the free surface. The nature of the structural changes to be anticipated depends ultimately 
upon the numerical value of a certain “internal” Froude Number which is a function of the 


initial velocity distribution as well as of the total density range. 


A. Introductory comments 


The immediate purpose of the theoretical 
analyses presented below is to suggest a me- 
chanism through which the momentum of 
certain stratified air and ocean currents may 
be concentrated vertically into a shallow 
layer or jet, and to establish a criterion from 
which the limiting vertical distribution of 
velocity in such a jet may be determined. 
Recent aerological investigations into the 
structure of the upper atmosphere over North 
America have shown that the prevailing upper 
westerlies of our latitudes, which may be 
said to constitute a “statistical” or climato- 
logical jet, frequently consist of two or three 


narrow zones of extremely high winds, 
stretched out along the lines of flow and se- 
parated by regions of weaker winds. It is with 
these individual jets, or streaks of high wind, 
that we wish to concern ourselves. One major 
aim of the investigation into the general circu- 
lation of the atmosphere of which the present 
analyses form a part, is to provide a satis- 
factory answer for this “streakiness”” which is 
becoming increasingly apparent as detailed 
synoptic data for the upper atmosphere 
accumulate. 

Through the courtesy of the Woods Hole 
Oceanographic Institution the writer has had 
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the opportunity to examine some of the data 
collected during a recent, broadly planned 
synoptic survey of the Gulf Stream system 
east of Cape Hatteras. Some of the results of 
this survey are described by FUGLISTER and 
WORTHINGTON (1951) in this issue of Tellus. 
These observations show that the current, 
after leaving the continental shelf of North 
America, and in spite of considerable entrain- 
ment of “slope” water along its left edge, 
remains narrow and well-defined, with a 
surface maximum speed of about five knots 
and a total width of approximately fifty 
nautical miles. In general, the behavior of 
the current during the survey period was 
strikingly similar to that of an atmospheric 
jet stream. Waves and meander-patterns were 
observed, including a complete cut-off of a 
cyclonic vortex on the Bermuda side of the 
stream. It is true that the length scale of these 
formations was considerably smaller than that 
of the corresponding phenomena in the at- 
mosphere. However, if it may be assumed 
that the horizontal meander-patterns are con- 
trolled by the planetary vorticity effect such 
‘a reduction in scale must be expected and the 


1 

ratio should have the value (ze) la where 
jet 
u, is the maximum speed of the Gulf Stream 
and wj+ the maximum speed of a fully de- 
veloped atmospheric jet. The observed ratio 
of length scales seems to be in fair accord 
with this theoretical conclusion. Hence, in 
view of the similarity between these apparently 
unrelated current systems it is not unreasonable 
to assume, as a working hypothesis, that the 
factors controlling the shape and behavior of 
jets must be fairly independent of their driving 
mechanism, and derivable from quite general 
dynamic principles. 

Recent theoretical investigations by Srom- 
MEL (1948), HiDAKaA (1949) and Munk (1950, 
1950 b) have shown that the broad-scale 
statistical features of the horizontal current 
distribution in the North Atlantic Ocean 
may be interpreted as the result of the 
permanent North Atlantic wind system. 
These highly successful analyses are based 
on the assumption that the vorticity of 
vertical water columns may be expressed as 
a balance between the changes in vorticity 
brought about through wind stresses, through 


latitude displacements, and through lateral 
mixing. The theory accounts in a remarkable 
manner for the large-scale features of the 
ocean circulation, in particular for the west- 
ward displacement of the North Atlantic 
eddy, but is incapable of accounting for such 
features as the presence of an extremely 
narrow meandering current in the western 
part of the ocean. 

This point of view is in accord with the 
results of an analysis undertaken by Mr 
H. Stommel of the Woods Hole Oceano- 
graphic Institution and presented as a sep- 
arate note in the present issue of Tellus 
(STOMMEL, 1951). Mr Stommel utilizes a 
composite map which contains the individual 
synoptic Gulf Stream patterns observed on 
various occasions during the survey referred to 
above, as well as individual positions of the 
current determined during earlier cruises 
(fig. 2 in FUGLISTER and WORTHINGTON, 1951).\ 
He then computes a mean value of the hori- 
zontal eddy viscosity (K) that would cor- 
respond to this statistical picture of the current 
distribution. The calculation -indicates that 
if the meander patterns are interpreted as a 
manifestation of lateral turbulence the asso- 
ciated eddy viscosity must have a value of 
about 2.108 cm? sec-1; thus it is about, five 
times as large as the value used by Munk in his 
theoretical analyses. This may seem to be a 
large discrepancy. However, the irregular and 
transient horizontal currents of the North 
Atlantic are certainly more intense in the Gulf 
Stream region than elsewhere. Furthermore, 
the horizontal eddy viscosity enters into 
Munk’s calculations in the form of a “wave” 
number which contains the third root of K. 
Hence the computed current patterns are 
relatively insensitive to inaccuracies in the as- 
sumed value of K. It may, therefore, be as- 


serted that Mr Stommel’s calculations lend — 


support to the statistical interpretation pro- 
posed above for the phenomena described by 
the Munk-Hidaka-Stommel theory. 
Similarly, the initial studies of the atmos- 
pheric jet stream (RossBy 1947, UNIVERSITY 
OF CHICAGO 1947) did not distinguish clearly 
between the permanent, climatological fea- 
tures of the middle latitude west wind belt on 
the one hand, and the details of transient local 
jets of extreme sharpness on the other hand. 
These first analyses were based on the as- 
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sumption that the basic features of the per- 
manent jet could be accounted for by quasi- 
horizontal exchange processes.! The meteoro- 
logical theories which may be outlined on the 
basis of this hypothesis account in a fair man- 
ner for the average strength of the middle 
latitude westerlies, and for the average dis- 
tribution of horizontal wind shear at the 
tropopause level, from polar latitudes south- 
ward to the vicinity of the center of the west 
wind belt. However, the rapid decrease of 
wind speed observed on the south side of 
the west wind maximum in synoptic cross 
sections through individual jets can hardly be 
accounted for by quasi-horizontal motions, 
and such features as the occurrence of highly 
concentrated double or multiple jets are not 
likely to be explained in terms of quasi-baro- 
tropic models. 


Such evidence as is now available indicates 
that strong individual atmospheric jets are 
limited both in length and in life span. Likewise 
there are reasons to believe that the concen- 
trated current found in the Woods Hole 
survey of the Gulf Stream may be limited 
in extent, and breaks up into several branches 
as the current reaches the Grand Banks. One 
is thus led to search for mechanisms capable 
of bringing about local concentrations of 
momentum within the broad “climatological” 
circulation patterns of the oceans and the 


1 To the extent that these analyses have made use of 
vorticity considerations they differ from the corre- 
sponding oceanographic theory in one basic aspect. 
Munk’s calculations consider the effect of lateral exchange 
processes on the distribution of the vertical component 
of the relative vorticity while the meteorological analyses 
assume that it is the vertical component of the absolute 
vorticity which is being redistributed by lateral mixing. 
Apparently because of the fact that the oceanographic 
problems considered until now have dealt with closed 
basins in low or middle latitudes the omission of the 
earth’s own vorticity component from the diffusion 
terms does not seem to invalidate the results of the cal- 
culations. . 

The kinetic energy and the vorticity of the large-scale 
quasi-horizontal eddy motions in the atmosphere must 
be derived from barocline processes; the distribution 
and intensity of the zonal motions established through 
such lateral exchange processes should be fairly indepen- 
dent of the specific nature of the mechanisms through 
which the lateral ‘“eddy’’ motion is created. In the ocean, 
the lateral exchange motion may derive some of its 
energy from the effect of stresses exerted by the non- 
permanent wind systems, but may also to a considerable 
extent be fed by the manifest meandering instability of 
such narrow currents as the Gulf Stream. 
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atmosphere. These local jets always appear to 
be associated with a strong solenoid concentra- 
tion, as required by the thermal wind equa- 
tion (circulation theorem) and with intense 
vertical shear in the fluid immediately below 
the jet, suggesting that the vertical stability 
distribution may plan an important role in 
the concentration of momentum. 


It should be stressed that the small width 
and the high speed of the upper layers of the 
Gulf Stream rather definitely exclude the 
possibility of interpreting the current in this 
part of the Atlantic as the result of momentum 
added locally by the prevailing winds. In this 
connection it is significant that the California 
current, which is considered as driven by 
the local wind systems off the Pacific Coast 
of North America is broader, more diffuse 
and far slower than the Gulf Stream. This 
comparison strengthens the impression that 
individual sharp jets result from a concentrating 
mechanism which is a basic characteristic of strati- 


fied air and ocean currents. 


A clue to the possible nature of the vertical 
concentrating mechanism may be contained 
in the results of certain aerological case studies 
which are being undertaken at the present time 
in Chicago. 

On several occasions during the summer 
of 1950 it was observed that westerly currents 
impinging on the Pacific Coast of North 
America and characterized by a fairly flat 
wind speed maximum near the 200-mb level, 
in crossing the mountains experienced à 
marked change in vertical structure. East 
of the mountains these westerly currents 
showed very slight air motion below the 
300-mb level. Above 300 mb, however, a 
sharp jet appeared, concentrated both ver- 
tically and horizontally into a maximum of 
much greater intensity than the one observed 
on the West Coast. It is hard to avoid the 
impression that the concentration must have 
been associated with the momentum losses in 
the lower portion of the air stream which 
resulted from the passage over the moun- 
tainous parts of the continent. 

On the basis of such clues it appears reason- 
able to attempt to ascertain what the vertical 
profile of velocity would be in an immiscible, 
stratified current which is subjected to mo- 
mentum losses through contact with the 
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underlying (rigid or fluid) boundary, or at 
its lateral boundaries, but required to trans- 
port, within each material fluid sheet, a pre- 
scribed amount of fluid per unit time. It 
would be particularly valuable to know the 
nature of the particular velocity profile which 
corresponds to a minimum value of the transfer 
of momentum per time across a vertical 
plane normal to the current axis. It is reason- 
able to assume that this profile represents a 
limiting state which would be approached 
gradually by any stratified current subjected 
to momentum losses but prevented from 
escaping to the sides. 

The withdrawal of momentum from the 
interior of a stratified current demands fric- 
tional (turbulent) interaction between layers. 
If we nevertheless are permitted to apply the 
principle of mass continuity to the individual 
fluid layers, this must imply that the individual 
fluid sheets remain immiscible, or that the 
eddy diffusion of mass properties is small in 
comparison with the eddy diffusion of mo- 
mentum. This is in good agreement with 
oceanographic experience, which indicates 
that in stable water masses the ratio of eddy 
conductivity or diffusivity to eddy viscosity 
may sink to o.1 or less. This point has been 
discussed by GOLDSTEIN (1938, pp. 229—231). 

The velocity profiles obtained from the 
principle outlined above are characterized by 
a rapid increase of current speed with in- 
creasing height above the underlying bound- 
ary, and also by a marked increase of vertical 
stability upward. When applied to well 
developed atmospheric jet streams the theory 
appears to lead to a reasonable value for the 
wind speed at the upper inversion. It may 
likewise be applied to well-defined ocean 
currents in which case it leads to surface 
current speeds in good agreement with the 
observational data collected during the recent 
Gulf Stream survey. 

The pattern for a theoretical determination 
of the profile of minimum momentum transfer 
is well known from fluid mechanics. In the 
case of flow of water in an open flume the 
transfer of momentum is given by 


D D? 
MT = f (gu? +p)dz = 0 (up as =): 


(1) 


D being the depth and u the velocity, assumed 
to be independent of depth. The water pressure 
p is computed from the hydrostatic equation 


and counted from the free surface. The 
volume transport V is given by 
V = 4D = constant (2) 


and hence 


pa D2 
ut =o (5+) () 


Variation shows that 
SMD—e(ep— 5) oD, ( 


and thus the extreme (minimum) value of 
MT occurs when 


Edie, EEE (5) 


If the flow takes place in a shallow layer D, 
flowing on top of a deep non-moving layer 
of density @ (> @), the corresponding value 
of the critical current velocity is given by 


0b = O O0 
u? = g — D = yD € = g eZe), 
Ob b 


(6) 


We shall now attempt to use the principle 
outlined above in an analysis of the vertical 
concentration of momentum in straight water 
currents in which density varies with depth. 
In part II, to be published later, we shall 
consider the vertical distribution of momentum 
in air currents and finally the mechanism 
which appears to be responsible for the hori- 
zontal concentration of momentum in air 
and ocean. 


B. Vertical concentration of momentum in 
water currents with a free surface 


We consider a straight, parallel stream 
incapable of escape motion sideways. The 
water is considered incompressible and the 
density varies with depth. Under these assump- 
tions, the transfer of momentum across a 
vertical strip normal to the current axis is 
given by 
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MT = 4 (ou? + p) dz, (7) 


where the vertical coordinate z is counted 
upward from the bottom and where u is 
current speed, p the water pressure (assumed 
to be hydrostatic). We shall now introduce 
a function of the density as the independent 
vertical coordinate, to help us incorporate 
and utilize the fact that the density is constant 
individually, 


d 
SEN (g 
£ 
We write 
@ = 9 (1 — 2x0) (9) 


where 04 is the density of the bottom water, 
and 


(10) 


is the surface density. Hence o, our new 
independent variable, ranges from o in the 
bottom layer to I at the sea surface, and 


Qs = % (T— 2%) 


= 0b — Os 
Ob 


2H 


(11) 


is the percentage range of density over the 
entire water column. 
In terms of the new variable, 


(12) 


I 
MT = f (ou? + p) 2 do, 
Oo 
where 


ALT (13) 


We shall next assume that the mass transport 
in every infinitesimal isopycnic layer remains 
constant during the variation process, though 
it may vary from layer to layer. Hence 


(14) 


where the subscript o refers to the initial con- 
ditions. This expression v for the mass trans- 
port may be used to eliminate u from MT 
and one finds 


ouzdo = ou, Z do = v(a)do, 


Because of the hydrostatic balance 


(16) 


| ioe x 
and thus one finally obtains 


I 


Mie [& u do. (17) 


R 08 


Our variation problem consists in determining 
the particular function p of o which reduces 
MT to a minimum for the prescribed distribu- 
tion of v with o. The variation of p vanishes at 
the sea surface and we shall assume that it 
vanishes also at great depths.? Under these 
circumstances the minimum value of MT is 
characterized by the fact that 


ö(MT) = 


a UND er ] “i. 
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(18) 


This is true for arbitrary values of öp provided 
the function p satisfies Euler’s equation 


Po E 2] = 
g dolp gl ” (19) 
which, after substitutions, reduces to 
Da “lle 
PAT oy 


Introducing the specific volume « instead of 
the density, one finds 


(20) 


| du? = pda | 


| (1) 


2 If the latter assumption is dropped a critical velocity 
corresponding to the speed of external gravity waves is 
obtained. This type of flow is clearly inconsistent with 
the prescribed external conditions of motion and it 
may be assumed that such variations in the bottom 
pressure quickly are dispersed through gravity waves. 
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To determine the final velocity distribution 
from the initial mass transport distribution it 
is necessary to combine (21) with (14), or 


(22) 


We may test the character of our basic 
result (21) by applying it to a simple, double 
current system in which a lighter fluid of 
depth D, and specific volume &, flows over 
a deep, immobile fluid of specific volume 
& (< &,). In this case du? and dx vanish every- 
where except at the internal boundary where 


oudz = v (o) do 


(23) 


2 


I 
Pr = 028 De (a = =) 
Hence, (21) reduces to 


Hy” = 02gDz (%,— 1) 


or 


in agreement with the result stated in (6). 

The basic equation (21) is particularly 
simple to evaluate for ocean currents, in 
view of the exceedingly small percentage 
variation of 0, which permits us to replace the 
variable o with a suitable mean value o in 
the computation of p. One then finds 


p = 082, 


where the vertical coordinate now is counted 
downward. 

Before attempting to determine theoretically 
examples of current profiles of minimum 
momentum transfer which correspond to 
prescribed initial volume transports we shall 
rewrite our basic equation so as to permit us 
to calculate the limiting velocity profiles 
which correspond to the observed density 
distribution in the central part of a fully 
developed ocean current. We have, with 
good approximation, 


du? do, 
plicated te io ==. 
ls d 


dz (26) 


where z increases downward and 6, repre- 
sents the quantity 


(27) 


In this expression 035,0,0 is the standard 
density corresponding to a salinity of 35 °/oo, 
atmospheric pressure and a temperature of 0° C. 
With sufficient accuracy (26) may be written 
in the form 


dy? _ do, 
Be (ead 


Gr = 1,000 (0 — 035,0, o)- 


(28) 


where z and u are in cm and cm sec! re- 
spectively. 

We shall at first attempt to represent the 
o,-distribution curve by the following ana- 
lytical expression: 


(29) 


It is evident that /\ measures the total range 
in o,. In the Gulf Stream region an adequate 
value for /\ appears to be 4.5. To determine w, 
note that (29) upon differentiation gives 


do, = 2.1] N» >» etz? 
dee MAL © 


O1 = Où — /\e 


(30) 


indicating that the o,-gradient disappears at 
the sea surface (within the homogeneous layer) 
and at great depths. The maximum gradient 


occurs at 


ats aa 


VY 2u 


(31) 


N 


which may be considered as the depth of the 
center of the thermocline. 
One then finds, from (28) 


W—u2=—2uA f ae dz, (32) 
[6] 


us being the current speed at sea level. Since u 
must vanish at great depths, 


ue = 2u /\f ze" dz = 
O 
Ix 
= /ND Va ED! (33) 
For A = 4.5 this gives 
He SOD SUD} 


(34) 
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Expressing u, in m.p.s. and D in meters 
this is equivalent to 


= 0,237 VD, 


(35) 
or, for a thermocline depth D of 100 m, 
4; == 2.37 mps, 


in good agreement with the observed maxi- 
mum surface speed of the Gulf Stream. For a 
D-value of 200 m the value of u, is raised to 
about 3.3 mps, somewhat in excess of the 
observed maximum surface speed. 

The vertical distribution of # obtained from 
(32) is such that the maximum vertical shear 
of the current occurs well below the center of 
the thermocline. 

The distribution of o, assumed above gives 
a more gradual density variation through the 
thermocline than normally observed in the 
Gulf Stream. We shall attempt to investigate 
the effect of a steep density gradient by 
assuming the same total range A in a, as 
before, but limited to a layer of thickness D. 
The thickness of the superimposed homo- 
geneous layer is H. Then 


do, 
7 =0 (36) 


for z < Hand forz > H + D, and 


re (37) 
ae. 


for H < z < H + D. It follows that 


ut =A(H+ =D). (38) 
Using 100 m for the depth of the homo- 
geneous layer and 200 m as the total thickness 
of the total thermocline one finds, for the 
same /\ as before, 


us = 3.0 mps. 


We shall next consider an initial value pro- 
blem. A uniformly stratified current of speed 
u, and depth D, flows on top of a homo- 
geneous bottom layer of density @, in which 


the volume transport vanishes, and is permitted , 


to readjust itself to a minimum momentum 


transfer current profile. We shall attempt to 
determine the distributions of density and 
velocity in the final state. Introducing the 
new variable 


D (1 + 2% G) (39) 
where 0, is the surface density and 


(40) 


the deep water density, it follows from the 
assumption of a uniform initial stratification 
that initially (subscript 0) 


z dz 
ni 


Hence the continuity requirement takes the 
form 


Op = @; (I F 2%) 


(41) 


Ua 


ee 
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(42) 


Our basic equation gives, with good ap- 
proximation, 


du? do 
Ae == ve FBS (y == 2 #9) (43) 
re do . , 
or, through elimination of 1 and integration, 
Et D? — 2°). (44) 


D D: 


Substitution in (42) gives 
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It follows from our basic equation that ¢ must 
be continuous wherever u is continuous. 
Hence, at the depth D, where u vanishes, o 


(45) 


‘must reach its maximum value of unity, 


corresponding to the bottom water density. 
Thus 


ur 
== S=S— —_— . gf 6 
Pen i in De 1 (4 ) 
and D may then be determined from 
Dia oft ig 
D; i re (47) 
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Fig. 1. Transformation of uniform current with constant 
vertical density gradient into flow characterized by a 
minimum value of the momentum transfer. The initial 
uniform velocity distribution is given by the heavy 
broken line, the final velocity distribution by the heavy 
full line. The density distributions before and after 
adjustment are given by lines marked ojnitial and Ofinal- 
Note that depth of final current is one half of initial 
depth. The total percentage density range has the value 
0.0045. 


It follows that the current must become 
shallower, and the bottom layer increase in 
thickness, whenever ty falls below the critical 
value #6 criticaa defined by 


Ug — Uo critical = Ve + > (48) 


The final surface current speed u may be 
computed from 1, by means of the expression 


3mps 


Ba Ds 
“L(Y (49) 


The process of adjustment is illustrated 
graphically in fig. 1 for A = 4.5, Do = 600 
m, #9 = 0.75 mps. In this case D = 300 m and 
== 2.25 pss 

It is clear that the non-dimensional quantity 


F, defined by 


og 
F= 5 D” (50) 


has the structure of a Froude number, in 
which the acceleration of gravity has been 
reduced in proportion to the total percentage 
density range of the Auid. The results obtained 
show that this new number determines the 
nature of the internal, or barocline, motions 
of the current, whereas the standard Froude 
number controls the external, or barotropic 
motions, such as long gravity waves at the 
free surface. 

The analysis presented above indicates that 
a uniformly stratified current, subjected to 
momentum losses through frictional action, 
will be concentrated to the lighter surface 
layers whenever the internal Froude number 
is less than a certain critical value (in this 
case 1/6). For Froude numbers in excess of 
this critical value momentum losses must 
result in an increase in depth. Also in this 
case some shearing motion develops, but it 
will be less than in the sub-critical case. 

It should be noted that whenever the final 
surface velocity u, exceeds the initial speed 
Uo, the free surface must slope downward in 
the direction of the current axis. The minimum 
slope required is easily determined from the 
Bernoulli equation. 

Though it is difficult to make exact deter- 
minations of the internal Froude number for 
real ocean currents with their arbitrary velocity 
profiles it would nevertheless appear, from 
rough estimates, that these currents generally 
are characterized by sub-critical values of F. 
Hence, when subjected to decelerations (mo- 
mentum losses), they would have a tendency 
to develop strong shearing motion with 
increasing velocities and increasing stability 
near the surface, decreasing velocities and 
stabilities further down. As far as the Gulf 
Stream is concerned, this would not necessarily 
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be the case in the Straits of Florida or over the 
Continental Shelf, where the range in O7 is 
smaller than further downstream and where 
no homogeneous deep water is available to 
facilitate separation of the current from the 
bottom. After the current leaves Cape Hatteras, 
however, the momentum it gains through 
direct wind action on the narrow strip exposed 
to the atmosphere is presumably incapable of 
balancing the losses which result from in- 
teraction with the deeper water masses and 
through lateral mixing. Hence the current 
tends to become more and more superficial; 
through this process the high surface velocities 
are maintained. 

The stagnant water masses which surround 
the swift moving upper layers of the current 
possess, as a consequence of the earth’s rota- 
tion, a high degree of inertia stability. They 
will therefore tend to resist changes in width 
of the current, in proportion to its superficial 
concentration. It follows that the vertical 
concentration of momentum must be accom- 
panied by a corresponding sharpening of the 
horizontal current profile. This effect will be 
discussed in Part II. There is, however, another 
possible consequence of the earth’s rotation 
which deserves mention at the present time, 
namely the effect on the current direction. 
As the upper layers shrink vertically and stretch 
out in the direction of motion they will 
necessarily acquire increased anticyclonic vor- 
ticity. They will therefore bend more and 
more sharply clockwise, away from shore, 
while the deeper strata, which are subjected 
to vertical stretching, tend to bend cyclonically 
toward the coast. It is reasonable to assume that 
this process must facilitate upwelling of colder 
water along the left edge of the current. The 
anticyclonic deflection may ultimately turn 
the current around a full right angle. This 
turning, in combination with the tendency 
toward upwelling along the left edge, may at 
least in part contribute to the formation of 
the peculiar, long-necked north-south meander 
pattern observed during the recent Gulf 
Stream survey. Similar patterns have been 
observed in the atmosphere, at the tropopause 
level, east of the Rocky Mountains. 

Our basic result up to this point, equation 
(21), is derived under the assumption that the 
mass transport does not vanish anywhere 
within the medium; the equation shows 


that the highest current speed then necessarily 
must be associated with the largest specific 
volumes. Hence the swiftest current must 
flow right at the surface. This would be 
true even in a current system which before 
adjustment had been characterized by a 
velocity maximum in the interior of the 
fluid. If the homogeneous top-layer before 
adjustment is characterized by a small but 
non-vanishing volume transport its thickness 
would decrease until the appropriate maximum 
speed. is reached. 

This result may be in accord with normally 
prevailing conditions, even though certain 
cross sections calculated in earlier years from 
the more sparse data then available show 
current maxima well below the free surface. 
This is the case, for instance, with the current 
values published by Iselin (1936) for a section 
across the Gulf Stream off Chesapeake Bay, 
in which a maximum of 120 cmps is located 
at a depth of about 200 m. New and detailed 
measurements are required to settle this point. 

The problem of current maxima below the 
free surface may conceivably be of some in- 
terest in connection with the analysis of river 
flow into harbors or tidal basins, where water 
masses of different density are brought into 
contact with each other. It is certainly of 
considerable meteorological importance. We 
shall therfore attempt to extend our theoretical 
analysis to include certain types of flow in 
which a stable maximum may develop also 
in the interior of the medium. 

It will be assumed that the moving, stratified 
fluid is overlain by a homogeneous layer, in 
which the volume transport vanishes. Un- 
derneath the current there is another homo- 
geneous layer of vanishing volume transport. 

Furthermore we shall assume that the specific 
volume & in the stratified current varies 
continuously from &,, the value in the homo- 
geneous bottom layer, to &, the value in the 
surface homogeneous layer. We introduce, 
as before, 


œ = ay (I +20), & = % (i +2), (51) 


and hence o runs from o at the bottom of the 
stratified layer to I at its top. 

In the analysis of the momentum transport 
of this stratified current, it can no longer be 
assumed that the pressure variations along 
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the upper and lower boundaries (a and b) 
disappear. Because of the absence of motion 
in the surroundings (apart from transient 
motions during adjustment), it follows that 


(52) 


It should be noted that under the assumptions 
made above, the level of the free surface must 
remain constant. This is readily seen from the 
fact that in the absence of all other motions 
but transient gravity waves no dissipation of 
energy can occur within the homogeneous 
top layer. Hence the Bernoulli equation may 
be applied to this layer; since the kinetic 
energy vanishes and the sea level pressure 
remains constant it follows that the height of 
the free surface can not change. For similar 
reasons it follows that the pressure at any 
fixed level within the lower homogeneous 
layer cannot change. 

To calculate the momentum transfer we 
extend our integral from some fixed level H, in 
the lower layer through the stratified layer up 
to another fixed, but otherwise arbitrary level 
H, in the upper homogeneous layer. Hence 
the variable portion of the momentum 
transport is contained in the expression 


Opa = — L0, 04, Opy = — gov Ob. 


H H 


a b a 
MT = / (ov? + p)dz =f pdz + [ pdz + 
ae ; dz 
ar J (ou? + p) zdo, (: = z) (53) 
or 
+ | — H, = 
MT = pu, (b— Hi) — ong = | 
Hi — a) 


ae fr (ou? + p) zdo. (54) 


It will next be shown that in computing 
the variation of MT from the expression given 
above one must take into account a restrictive 
condition associated with the lack of motion 
in the two homogeneous layers. 

There can be no energy losses in these 
two media since they remain motionless. Since 
they, furthermore, are incompressible and 


~ 


homogeneous it follows that the energy 
per unit mass (E) is constant from point to 
point and with time in each of these two 
strata. One finds 


(55) 
(56) 


where h is the constant height of the free 
surface above the bottom (z = 0) and poo the 
constant bottom pressure. In the intermediate 
stratified current E varies with height and has 
the value 


E, = px, + gz = gh 
and 
Ey. = pay + gz = % Por 


Be 


2 (57) 


= PO Gee 


The change from layer to layer is given by 


dE 1 du? do dp dz 


ES PAL DU ducs 


The two last terms on the right side of this 
equation vanish in consequence of the hydro- 
static requirement. Hence 


dE 7 du? dx 


do 2 do se: 


(59) 


and, if this equation be integrated across the 
current, 


— gh — & Poo = Constant. 


(60) 


Following standard procedures we multiply 
the constant integral for E,— E, with an 
arbitrary parameter À and add to the expression 
for MT. The variation of the sum is easily 
computed. Keeping in mind that | 


Ph — PH, — HL (b Zu; Ah), Pare [ist + 
+ 08 (Ha — a), (61) 
one finds 
Opa Opp 
en + 
ges Pgo 


(62) 


6 LMT AE ESS) 
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+ 6S [ou2z + pz + 2+ 2x0, p] do. 


[e] 
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For extreme values of momentum transfer this 
gives, following the same procedure as before, 


o — PaO Pa _ 
8 Qs 


I 
Ss [+ 3-2) dp + 
Se op | de 
teem gor 1” 


where v represents the mass transport per unit 
interval of o. The integral may be evaluated 
by parts and one then finds 


0 2 ‘ a 
en [CE — 2) op | + 
8% RB 80 b 
x . 
| 2 3 Je gus 2 ate 
+ f [2x2 go 
o 


ENT. 
To = =) öpdo = o. 


Replacing v by the current speed one finds, 
after a few reductions, as condition for an 
extreme value of the momentum transfer 


Popo, 
analy 


(63) 


Pad Pa 


80s 


; [u2 dp], + 
d 
+ [ra -2) | öpdo — 0 
x (65) 
Hence, 
Ua ©, t= 0 (66) 
and 
a = (p + gd) 2x% (67) 
or | 
du? = (p + gd) da (68) 


Since u, and #, vanish, one may determine À by 
integration of this equation across the current. 


One finds 


er p 
= == da=— =, 
VER 5 


(69) 
£ en 


where p is the arithmetic mean value of the 
pressure with respect to specific volume. Hence 


du? = (p—p) da. (70) 


The modification so derived of our previous 
criterion is perhaps of greater potential use- 
fulness when applied to the atmosphere, for 
instance in the analysis of cold air outbreaks 
in the lower troposphere. Nevertheless; a 
simple application to incompressible media 
will be given. It will be assumed that initially, 
within the stratified layer, 


dz 
u = Uy = constant, — = 
oO 


(OE. 


(71) 


where D, is the initial thickness of the stratified 
layer. Hence, in the final state 


dz 
QU Te = 0m D; (72) 


and substitution of this expression in (70) gives 
(73) 


To integrate this equation we make use of 
the fact that to a very good degree of approxi- 
mation 


Pp) = g & —2), (74) 


co (p 


where Z is the level of p and the vertical 
coordinate z is counted upward. Hence 


du y (2 — 2) 


dz 24,0% 


(75) 


and after integration, 


22 — bh? 
ae zn E a 2 =|. (76) 


b as before being the lower boundary. # must 
vanish at the upper boundary a and hence 


Z (a—b) = Ma (77) 
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Of 


a +b 


=: (78) 


>= 


A critical Froude number may now be 
obtained from the requirement that o must 
vary from 0 to 1 as z goes from b to a. This 
may be seen from the calculation presented 
below. 

If one introduces a vertical coordinate y, 
measured upward from b, the expression for 
u given in (76) may be changed into 


= ue ( 


u— 4m D, U (79) 


D—y), 


where D is the final depth of the stratified 
current, 


ss D 
Den hi. 


(80) 


The expression for de follows from a combi- 


dy 
nation of (72) and (79). One finds 


do y 
dy 4? Dy 


39 (D —y), (81) 


and hence, after integration 


k Y y? D ] 
i: 4 Up? Do? [ 2 anh 22 


At the top of the stratified layer (y = D) this 
must have the value 1. Hence 


DA te Ash DA 
RATE Ts (83) 
and 
DNS Hee 
(D) = rb, a 


The critical Froude number now has the 
value 1/24. When 
tee et 


F= < 
7 Do 24 


(85) 


the stratified current will contract and speed 
up at its center. For F > 1/24 the current will 
expand its cross section and slow down, even 


though there still will be a concentration of 
motion at its center. One may introduce an 
average value of the stability for the stratified 
layer in its initial state, 


__ gl m) 2g _ y 
& &h ID). ID); ID, | 


Then 
Up 


HE gs Di? 


and the critical relationship becomes 


(88) 


[gs 
Uo == ID) \ DE 


Further discussions of the velocity profiles 
derived ‘above must be postponed until 
observational data are available. 

Professor A. Craya, of Grenoble, with 
whom the author has had the opportunity to 
discuss the basic principles of this paper, has 
suggested that the restrictive condition (60) 
derived above has a much broader validity 
than here indicated. Indeed, from the hydro- 
static condition it follows that 


gh— a poo = SI pdz, (89) 


ab 


where h is the height of the free surface 
above the bottom, po the constant bottom 
pressure and &, and a, the specific volumes 
at the bottom and the free surface respec- 
tively. Since the value of h is very nearly, 
even though not completely, determined by 
the depth of the surrounding flaid, it follows 
that the integral on the right side of (80) 
must be kept constant during the variation 
of MT. Professor Craya’s analysis is pres- 
ented in a separate contribution to this issue 
of Tellus (Craya, 1951). 

We shall finally rewrite our basic equation 
(21) through the introduction of the vertical 
stability s. Counting z positive upward one 


finds 


and thus, through substitution in (21), 
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JIS — En (91) 


It will be shown in part II of this study that 
in spite of the compressibility of air the same 
equation (91) is valid also in the atmosphere. 


du? M: 
0 
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Critical Regimes of Flows with Density Stratification 


By A. CRAYA!, Iowa Institute of Hydraulic Researh 


(Manuscript received 16 January 1951) 


Abstract 


This study first discusses the analogy between flow in an open channel and either internal 
or surface currents in large bodies of fluid and shows that the concept of the critical regime 
is generally applicable to such cases. The critical regime of gravity currents due to density 
stratification is then examined in terms of both the momentum principle and the somewhat 
simpler energy principle; contrary to the case of open-channel hydraulics, the two principles 
result in different flow criteria. In the case of surface currents it is shown that, for constant 
piezometric or energy level at the surface, the critical regime corresponds to a velocity 
distribution which reaches a maximum some distance from the surface. 


In the foregoing paper, “On the Vertical 
and Horizontal Concentration of Momentum 
in Air and Ocean Currents’, Dr Rossby 
makes novel and significant use of the 
momentum principle, which has long been 
used in its more elementary form as one of 
the basic relationships of fluid mechanics. 
The central theme of his analysis, in fact, is 
the generalization of the concept of “‘mini- 
mum” conditions of momentum flux which 
is fundamental to open-channel hydraulics. 
As a result, not only does it open promising 
new horizons for hydraulics, but this field in 
turn may eventually be able to contribute to 
the extension of Dr Rossby’s original analysis 
— particularly in the parallel application of 
the equally fundamental energy principle. The 
writer, as a hydraulician, has hence sought to 
organize in a manner which is equally useful 
to each of the fields the pertinent aspects of 
both principles, and to present some supple- 
mentary results concerning surface currents. 

It is at once evident, of course, that the 
transition from the one degree of freedom of 
schematic open-channel flow to the many 
degrees of freedom of meteorology and 
oceanography introduces unavoidable difficul- 


t On leave from Etablissements Neyrpic and Uni- 
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ties and obscurities. I will therefore ask the 
reader’s permission to proceed very gradually 
and to recall briefly certain familiar fundamen- 
tal notions from the hydraulics of open chan- 
nels. The analogy between the flow in an open 
channel and the movement of density currents 
becomes more clear if we consider first the 
case of homogeneous gravity currents. We 
will then analyze the principal object of this 
study, that of the critical regime of stratified 
gravity currents. 


I. The critical regime in the 
hydraulics of open channels 


I. Hypothesis 


As it is well known, the concept of a critical 
regime is met essentially in the study of steady 
flow in open channels. Steadiness naturally 
implies the constancy of discharge from sec- 
tion to section. Henceforth we shall con- 
sider two-dimensional flow (infinitely large 
channel) where q shall denote the flow per 
unit width. 

An important characteristic of flow in 
open channels is the quasihydrostatic pressure 
distribution in a given transverse section, i.e., 
p + egz = constant. This assumption is based 
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on the hypothesis that the curvature of the 
streamlines is negligible. 

A further simplification employed in hy- 
draulics is that of a uniform velocity distribu- 
tion at a section, which thus is equal to the 
mean velocity of the flow 


q 
gern 
en 


2. Use of energy equation 


The counterpart of the energy principle for 
incompressible fluids is Bernoulli’s theorem, 
which states that (for a perfect fluid only) the 
quantity p + ogz + ou?/2 is constant along 
the same streamline. 

In view of the assumptions made, this 
quantity is the same for all the streamlines in 
a given section. IE we take the floor of the 
channel as the plane of reference, we obtain, 
after dividing by og, the total head 


q 
2 gh? 


This quantity varies from section to section of 
a channel, due to the changes of elevation of 
the floor and due to energy losses, but for a 
given discharge it cannot decrease below a 
minimum value which occurs for the critical 


depth defined by 
gh? = 
q? 


I 


In other words, a given discharge is allowed 
to low either at great depth and small velocity 
or at small depth and great velocity and the 
critical state corresponds to the intermediate 
stage when the depth equals twice the velocity 
head. For a channel with horizontal floor 
where the total head is constantly decreasing, 
there is then a tendency towards critical 
conditions. 


3. Use of the momentum principle 


Let us consider the integral across a section 
of the quantity 


p + ew 


The variation of this quantity from section to 
section represents, in view of the momentum 
principle, the horizontal component of the 
forces which the outer portions exert on the 
periphery of the flow between the two sec- 
tions. This at least is true if we neglect the 
normal forces of viscosity and the turbulent 
flux of momentum across the transverse 
sections. 

For the flow in a channel we obtain the 
quantity (which, for the sake of brevity, will 
be named henceforth the thrust) 


| he 2 
= E | ie) 


For a given discharge there is a minimum value 
of this thrust, which corresponds to the same 
critical depth defined previously: 


3 
RN 
q 


In a horizontal channel the thrust, according 
to its mechanical significance, cannot but 
decrease because of friction forces on the 
bottom, and we find again and for a slightly 
different reason a tendency toward critical 
conditions. 


4. The influence of unequal velocity 
distribution 


If we take into account the unequal velocity 
distribution, we have to express the mean total 


head as 


h 


CE) 2 
ehr [MER 


Ss 


[e) 


and the thrust as 


If, for example, similar velocity curves are 
considered, « and B will be independent of h 
and the critical regimes corresponding to the 
total head and the thrust will be, respectively, 
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Fig. 1. Gravity currents; (a) submerged, (b) surface. 


The two lines of attack, then, no longer lead 
to the same result and we can add that for 
the particular case under consideration B > x. 


II. Extension to homogenous gravity 
currents 


5. Preliminary remarks 


In what follows, we shall denote by “sub- 
merged gravity current” a flow confined be- 
tween two layers of fluid at rest with densities 
0, and @,, respectively (fig. 1 a). By “surface 
gravity current” shall be understood a flow 
with its free surface exposed to the atmosphere 
with pressure p = o and taking place on top 
of a fluid at rest with density 0, (fig. 1 b). 

The homogeneous gravity currents which 
we intend to examine first are those for which 
the layer in movement has a uniform density. 

For the submerged as well as the surface 
homogeneous gravity current, it can be seen 
from fig. 2 that the limits A and B of the 
‘current are not independent (the figures 
represent a hydrostatic pressure distribution, 
and the slope — og of AB is fixed). 

The interfaces A and B lie on either side 
of the fixed level Q and coincide with this 
level only when the current is reduced to an 
infinitely thin layer. The piezometric head 
line of cach layer is by definition the fictitious 
free surface, (p = 0) which could be obtained 


by prolonging this layer beyond its actual 
limits. Thus in fig. 2a the piezometric head 
line for the liquid at the bottom of density 
0 is the point No, for the intermediate liquid 
it is N, and for the surface layer it is N,. We 
can see from fig. 2a that the variation NN 
of the piezometric head line around the 
density discontinuity in B, for example, is 


equal to (po/g) A (1/0) or (4e/e) b No. The 
variations of the piezometric head line are 
then very small for very small variations of 
density, but it must be noted that they nev- 
ertheless command movements of great 
amplitude of the current interface. 


6. Minimum of energy flux 


As in open channels and with reference to 
the theorem of conservation of energy, we 
study the flux of the quantity; p + ogz + 
+ ou?/2 = E across a vertical section of the 
current. 

As E is a constant across the gravity current, 
its total value is given by gE. This flux con- 
tinuously decreases from section to section 
by the amount of the internal energy losses. 

(a) Let us consider first the case of the 
submerged flow. Fig. 2 indicates an imme- 
diate geometric relation between N,P (con- 
stant value of p + ogz across the current) and 
the thickness h of this current. We can express 
this relation by writing for the upper interface 


P + 092 = Pi 0021 = fi 7 Oz 
+ 921 (0 — 01) 
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and for the lower interface 


p + 082 = Po + 0820 = Po + 0092 + 


+ 92 @— 04) 


Taking p +087 = BZ, Po + 0082, = Zo; 
Pi + 01821 = 8Zı 
TE ioe 


A 26 = 
0 — 01 
~ = 


Z— Zo 


we obtain 
0. —@9 


from which we obtain 


90 — 01 


Z + constant 
0o— 0) (0— 0) 


er “=, = h = 
( 


Since E is defined (by the arbitrary datum 
plane) except for a constant, and since we are 
interested only in its variation, we can definitely 


adopt 


We find consequently for the critical regime 
which makes E a minimum 


(00 — 0) ae), ur 


2 


0 (00 — 01) qd 


Le ee | 
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. Pressure distribution in homogeneous currents; (a) submerged, (b) surface. 


; ; a 
and in the particular case of o = za 
Ao gh? 
0 


(b) The analysis of a surface current proceeds 
in an analogous manner with the help of 
figs D: 

We have as previously for the lower in- 
terface 


092 = Po + 0820 = Po + 00820 + 
+ 820 (0 — 00) 


pa 


which, by setting p + ogz = gZ and po + 
+ 00820 = gZo, yields 

Z— Lo 

0 — Co 


Es = 
< 


As far as the free surface is concerned, the 
relation is simply p + egz = 0gz,, which gives 
A == 20: 


From this we deduce 


Co + Z + constant 


Zo 2, = = 
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and the expression for E is then 


2 


19 


0 (00 — 0) 0 q 
E = = gh + 
00 8 2h 
which for the critical regime becomes, with 
Po 0 = 40, 

| degh 


— —— =] 


0o 49° 


It may be remarked that this criterion is the 
same as the one for the half part of a sub- 
merged symmetrical current. 


7. Minimum of thrust 


With reference to the momentum theorem, 
we shall analyze the integral of p + ou? 
across a vertical plane. This integral con- 
tinuously decreases from section to section 
by the amount of the longitudinal shear on 
the periphery of the current. 

If the lower liquid extends to infinity, the 
integral of the pressure will be infinite, but 
we can always subtract from it an invariable 
quantity which comprises the infinite portion. 

(a) For a submerged flow (fig. 2a), it is 
natural to subtract the constant distribution 
of pressure N,AQBR, so that for the variable 
part only the area AQB remains, which has 
the value 


OS hy + hy 
But 
= gh, £ho_ # 
23 I K T 


ax (hi ar hy) (go — @) (e — 0) 


00 — O1 


so that we obtain for the thrust — 


p = (@o—e) (e— a) : h? Le q° 
00 — O1 a) h 


and we get for the critical regime the same 
result as before. 

(b) In the same manner, for a surface current 
we shall subtract the fixed pressure distribu- 


tion OBR (fig. 2 b), so that the variable part 
of the pressure integral shall be reduced to 
the area AOB, i.e., 


But 


and the thrust can be written as 


e(o—eo) fh gq 
2 Os EUR TR 


and the agreement with the critical regime 
conditions obtained from energy considera- 
tions is again maintained. 


8. Remarks on surface currents 


We have already stressed the fact that in 
the case of homogeneous currents, A and B 
necessarily vary together (fig. 2b). With the 
very slight differences of density, and reason- 
able depths encountered in practice, the varia- 
tion of A is very small; but this reduced 
degree of freedom of A is the necessary 
counterpart of the far larger degree of freedom 
et à. 

Thus, if we rigidly fix the free surface A, 
there is but one regime possible. 

The same thing is true if we fix the total 
head of the topmost layer. In fact, on the 
basis of the assumption of uniform velocity 
distribution, we impose the condition of the 
constancy of the entire flow energy — i.e., 
according to previous calculations, 

2 
m = € 
Qo 2 gh? 
Thus there are at most 2 regimes corresponding 
to the given conditions. 

We shall see that things are different for 
stratified currents; this is due to the fact that 
they can vary in a more flexible and less 
monolithic manner than in the previous 
schemes. 
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Fig. 3. Pressure distribution in stratified currents; (a) submerged, (b) surface. 


III. Stratified gravity currents; critical 
regime of minimum energy 


9. Preliminary mathematical notions 


In the course of this analysis we shall have 
to consider the first variation of an integral 
of the form 


T= Jante 2.2752"). dx 


if Z is given a variation 6Z(x) in the fixed 
interval xox. 

In generalizing the more classical result 
pertaining to the case where the function 
contains only the first derivative Z’, we obtain 


easily 
a1 = | (35 — 2) 2]. + 
n E oz’ | an 
je + £a 
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If we seek a solution which will make I 
stationary for the conditions with the given 
limits, the solution must first of all satisfy the 
condition 


Op di ap Rd“ op 
IZ EST on? (1) 


The first variation is then reduced to 


dp d oo x1 
NS [GE =: sn) oz] " 
| OG) an] 
Ze Ee ÔZ |: (2) 


Depending upon the nature of the conditions 
imposed on Z and Z at the extremities of the 
interval, we can deduce, by cancelling this 
residual variation, the four constants of inte- 
gration of the differential Eq. (1). 


10. Condition of minimum of energy flux 


Generalizing the previous considerations, 
we are lead to examine for a stratified current 
the energy flux across a section: 
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As the density is constant in each thin layer 
fig. 1a), as well as the elementary discharge 
dq of the layer, we can set 


dg=udz=—k(o)do 6) 


The function k (0), which imposes a fixed 
distribution of the discharge between different 
layers, generalizes thus the notion of total 
given. discharge in the preceding calculations. 

We shall take as a characteristic distribution 
function of pressures, densities, and velocities 
in a section the quantity: 


gZ =p + ogz (4) 


This expression is closely related to the notion 
of piezometric level and like it, varies rela- 
tively little but commands the whole internal 
structure of the current. 

By differentiating Z we obtain 


gdZ = dp + ogdz + gzdo 


As dp + og dz = 0 from the assumption of 
hydrostatic distribution, we have 


dz 7 (s) 


de 


Consequently, by virtue of Eq. (4) the pressure 
is given by 


p =g(Z— eZ’ 


On the other hand, because of Eqs. (3) and 
(5), the velocity is given by 


The energy flux is written then in final form 
with this variable Z 


09 


k3 
le Harz + = za) do 


Q1 


and in order that this flux becomes stationary 
(in our case a minimum) Z must satisfy the 
condition 


HE ORNE 
t— i (Ba) = 


Te 
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ipa (218) = — ek () () 


II. Profiles of critical velocity 


In order to introduce dimensionless variables 
we shall set 


0 — O1 


and k (9) = 1/0 


where qo designates the total discharge and A0 
the difference of density 09 — 0, between the 
extreme layers of the current. 

Then, the differential equation (6) of the 
critical state becomes 


dè 
a (er) = —gaAof (r) 


Introducing also two successive integrals of f 


LE T 


A) So dr ® (r) = S F(r) dr 


ie) 


and a reference velocity w = V (40/0) ¢ qo 
(which is practically a constant if density 
varies very little) we have 


The function f (r) being positive, ® (r) is con- 


cave upwards and from fig. 4 we get a general 
picture of the required velocity distributions 
for critical conditions. Depending on the 
values of & and ß, the velocity curve con- 
tinuously increases from bottom to top, or 
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Fig. 4. Critical distributions of velocity. 


on the contrary has an intermediate maximum. 
It remains to eliminate this arbitrariness 
according to the conditions at the limits. 


12. The case of a submerged current 


The variation ö] for distribution of critical 
velocity can be written according to Eq. (2) 
and the value of 99/9Z" = — ok/Z"3 = ou? 
in the form 


20 


d 80 
108 rs 3 yo 4 
öl = le LAN + [es OZ le 


For a submerged current Z, and Z, (or the 
extreme piezometric levels) are imposed; 
consequently 6Z, and ÔZ; are zero, whereas 
values of 6Z’ = Ôz at both extremities are 
free and independent, (i.c., there is no rela- 
tion here between the displacement of A and B 
as in the case for homogeneous current). Then 
we must have u, = o and u, = 0, and in fig. 4, 
ar + ß is represented by the straight line D, 
which passes through A and B. 

The constants & and ß are then 


and 


The distribution of densities as a function of 
the height results immediately in 


w | Br 
qo Var + B— G(r) 


The total height h of the current is thus given 
in a general manner after setting JT = (Ao/e) 


(g h/q?): 


a ie: 
mg ee 


and in particular for a submerged current is 
given by 


yee er a 
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For an initial state of constant velocity and 
linear density distribution we have f = 1, 


Or) 7? /2,cand 
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This last value of the parameter JJ is calculated 
either by means of elliptic functions or by 
numerical integration. The value of JZ thus 
found is ZI = 17. 

It is also valuable to derive on hand of this 
general theory the results for the homo- 
geneous currents obtained more directly at 
the beginning. If q is the discharge between 
the layer of density 0, at the upper interface 
and a layer of density 0, we have 


Let us consider the case where this discharge 
is zero from 0; to 04, becomes qo for pa, and 
remains qo from 040 0,. This can be considered 
to be the same as the case of a homogeneous 
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current where the layers comprised between 
o, and 0, and between pa and 0, are like two 
infinitely thin sheets constituting the surfaces 
of discontinuity of density. 

The function f(r) is then zero everywhere 
except for the value of ry corresponding to 
oa. The function © is represented by the broken 
line AIB of fig. 5 and the value of the para- 
meter becomes 


a I 2 
men M4 
I 
- V ra Ge) 
This gives 


(e — 01) (00 — 0) gh° 
e(Qo— 01) q 


which is in complete agreement with the 
previous result. 


13. Case of a surface current 


Consider again the variation I for a distri- 
bution of critical velocity: 


A E (oui) oz |” a [ ew 2 |" 
do Q1 


Q1 


For the case of a surface current, the piezo- 
metric level of the lower liquid at rest is always 
imposed. Thus, as previously, Z, is given 
and ÔZ, is zero. 

The condition of the surface; however, is 
different. In fact, the free surface is characte- 
terized by pı = 0; 


hence, gZ, = pi + 0182 = 018 21 


07: = 01021 


l'E, 


Furthermore, as shown before, Z’ = dZ/do = 
= z, so that the expression for the variation 
öI becomes 


OT = pou dZ% + 2 E (eur) | En hoz 
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Fig. 5. Velocity diagram for homogenous current. 


The conditions that determine « and ß be- 
come then for the present problem 


te = 01 and) 42 = E (our) | 
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When expressed in terms of the relative vari- 
ables previously introduced, this last relation 


becomes 
20 [ d fu | 
dr WwW 1 
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The first condition u, = o requires that 
a + B = @(t), from which we deduce that 


no rage 
Hal a 


It follows, then, that for a surface current the 
distribution of critical velocity at first in- 
creases with elevation above the lower inter- 
face, passes through a maximum below the 
free surface, and then decreases (see fig. 4, 
where «r + ß is represented this time by the 
straight line D,). 

It must be added, however, that this rigorous 
result (there are no approximations in the 
preceding calculations) reduces practically for 
very small differences of density to a maximum 
velocity at the free surface; the slope « is then 
in fact very small and we have essentially 


Bela) %=0 
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Under these conditions the critical depth 
satisfies 


HP ered 
2 lap 


In particular for the case of f = 1 and ® = r?/2 
we obtain 


E I 
Bl dr Care dr 
V IT — Fur —— ia [ = 
/ı r yi—r 
& ate 2 3 


After integrating we get 


I 
Tea 
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We could also determine, by exactly the same 
method, the conditions of critical regime for 
a homogeneous current established differently 
at the beginning of this paper. 


14. Surface currents — supplementary 


remarks 


For surface currents a question comes up 
which does not intervene in the study of 
submerged currents; in the latter case the 
physical condition of a current limited on 
either side by liquid at rest is perfectly clear 
and can be unequivocally translated into 
mathematical conditions. For a surface current 
the same thing is true for the lower interface, 
but the question can be raised whether all 
required conditions for the free surface were 
properly taken into account. 

In the preceding calculations the free surface 
was left perfectly free to move along the ver- 
tical. Is this always true and will not other 
physical conditions limit this freedom of 
movement? One could think, for example of 
marine currents like the Gulf Stream. The 
presence of a mass of water at rest on either 
side of the current of limited width would 
cause a lateral entrainment, if the level of the 
current was allowed to vary. On the other 


hand, even for a two-dimensional flow 
(laboratory experiment in wide channels), 
if the drag on the topmost sheet is to be zero, 
the total energy of this sheet 2, + m?/2g 
must remain constant. 

The problem then becomes, to paraphrase 
Kipling, another story, which would require 
further developments. We want to note here 
simply that in the case of stratified currents 
we have the possibility of taking into account 
such a condition, whereas for a homogeneous 
current a condition imposed on the free sur- 
face would determine completely the flow. 

To show this for a relatively simple case let 
us imagine, for example, that the imposed 
condition is a free horizontal surface. We have 
then z, given and öz, = o. The freedom which 
remains for dz requires for the minimum 
of the energy flux that #, = 0; i.e. (Eq. 7) 
x + B = ®(1) and hence 


(4) -2(—-n +9) 0) 


w 


It is then possible to determine the correspond- 
ence between « and z (a calculation which 
seems useless to develop here in detail). 
From fig. 4 we can see that now we have 
a straight line D (representing «r + B) passing 
always through B, but susceptible of taking 
different slopes « depending on the condition 
of level imposed. The corresponding distri- 
bution of velocity has this time a net maximum 
under the free surface. The study of the case 
where the energy of the topmost sheet is 
imposed would led to analogous conclusions. 


IV. Stratified gravity currents; critical 
regime of minimum thrust 


15. Variation of the thrust 


Changing now the angle of attack, we take 
up the principle of momentum and conse- 
quently the study of the thrust. As for the 
homogeneous currents, we are led to consider 
as the variable part of the pressure integral 


Q1 
the cross-hatched area of fig. 3, Jp dz —F 
Q 


where F designates the area aAQBb for a 
submerged current (fig. 3 a) and the area OBb 
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for a surface current (fig. 3 b). We shall 
write the thrust then as 


01 
P=/f (p + ou?) dz—F 
00 


Taking again as the principal distribution 
function the piezometric level, or the quantity 
related to it, 


Ba 
the expression for the thrust becomes 


[21 
2 
pe f{[s22"- ogZ'Z" +0 = do—F 
20 


The distribution Z which will make P a 
minimum must first satisty Eq. (1) of section 9; 
that is, after simplifying, 


d? k? # 
en Û 


For the functions satisfying this differential 
equation the variation of P is then (Eq. 2, 
section 9) 
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From fig. 3 a, it is clear that for a submerged 
current we have for öF 


OF = p, 02; — pod Zo 
For a surface current it will be (Fig. 3 b) 
OF = — pydZ 
On. the other hand, for a submerged current 
Z, and Z, are imposed, hence Z,= 0, Z,= 0, 


and we have 


OP = [ourdz je 


For a surface current pı = 0; Zı = 01921, and Zo 
is given; the remaining variation 1s then 


d 
OP = Qo OZ + 01 lee (ou?) — | 024 
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16. Condition of minimum thrust 


The first thing that commands attention 
here is evidently that considerations of the 
minimum energy flux and of the minimum 
thrust do not lead to the same result. This 
will seem less surprising, however, if we recall 
that in the hydraulics of open channels this 
agreement is obtained only for the case of 
uniform velocity distribution. Such is not 
the case for stratified density currents where 
the velocity profile varies in accordance with 
the law of conservation of flow for each sheet. 

We can give to Eq. (8) different interesting 
forms. 

(a) If we set ou? = Y, we have 


If o varies very little, it is quite legitimate to 
use the following approximation. 
du? k (o) 


I au na = — 
do? Sg 0 


Introducing relative variables by 


we obtain 


We recall that the corresponding equation for 

the energy flux was 
dè &3 
‚de — 


I 


The present result is thus more difficult to 
analyse for the general case. 
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(b) We can also integrate Eq. (8) once and 
write: 


qe (oH) = ge +2) 6) 


Assuming again that @ varies very little and 
using the continuity equation, we get 


e kl) Ÿ 


— 2 


= g(z + a) 


For the particular case when k = 4/10 
(current with an initial uniform velocity 
and constant density gradient) we have 


du I Ae 


= BD g(z +) (to) 
(c) Lastly, if we integrate Eq. (8) twice, we 
have the alternative form 


(11) 


Written in the form u? = p/o + gz + gx — 
—pm/o and differentiated this equation be- 
comes, after taking into account also the 
hydrostatic equilibrium, 


QU? = p + 082 + 08% — Pm 


dur = a 
ee 


which was used by Professor Rossby in this 


special form. 


17. Critical depths 


From now on we are going to limit ourselves 
to the particular case of k(o) = 4/40, which 
alone is solvable without excessive calculations. 

(a) For a submerged current the conditions 
at the limits are #, = 0, u, = 0. Relation (10) 
shows that the distribution of velocities is 
parabolic. By taking the middle of the current 
as the origin of the 2’s, this distribution can 
be written 


he 
PR AL. «(= =) 
qo 0 4 
and the total flow satisfies 


+ h/e 
J udz= 9 


—hl2 


From here it follows immediately that 
Ag hs 


IR = 
0 fe 


This origin of z's corresponds to « = 0. 

Now Eq. (11) imposes on the two limits of 
the current where the velocity is zero the 
conditions 


8Zı + LAC — Pm = 0 
gZo + A 00 —Pm = O 
from which 
Zi — Zo 


01 0 


Consequently, the middle of the current is 
the juncture of the two straight lines repre- 
senting the hydrostatic pressure distributions 
imposed by fluids at rest surrounding the den- 
sity current, i.e., point Q of fig. 3 a. 

(b) Let us consider now a surface current 
and let us take as origin of ordinates z the 
piezometric level of the motionless liquid 
below (fig. 3 b). 

From section 15 the actual conditions of 
minimum are 


Uo = 0 
do 


Equation (10) gives for the free surface (p = 0) 
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O1ty? = 01821 + 018% — Pin 


and for the lower interface (p + 0 gz = 


Coto” = 008% — Pm = O 
On the other hand, from Eq. (9) 


[5 ew] =a +4) 


We have thus pm = 0 and «-= o. The velocity 
distribution is then, with regard to Eq. (10), 


u = — — — g(z?— 25°) 
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We see that the maximum of velocities is 
not exactly on the free surface (z = 21), but a 
little below, on the piezometric level z = 0 
of the motionless fluid (point Q of the figure) ; 
this will be a negligible difference if the 
density varies very little. Consequently we 
shall take z = o for the level of the free 

O 
surface and, taking / udz = qo, we obtain 


Zo 


3 
md og 
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If now we compare the result obtained for 
the minimum of thrust and that obtained for 
the minimum of energy flux, we find that 
the second criterion gives smaller critical 
values of IT. 

At the beginning of this study we performed 
an analogous comparison for the hydraulics 
of open channels and obtained opposite results. 
It must be added, however, that the two 
cases are not comparable, among other 
things because in the case of open channels 
we had assumed that the velocity curves 
remained similar during the variations. 


18. Supplementary remarks on surface 


currents 


One can also ask the same question as for the 
minimum of flux energy: Do we have to let 
the vertical displacements of the free surface 
be completely arbitrary during the variations 
or shall we subject them to some such con- 
ditions as the constancy of total head, 


Uy" 
z,+—=d (12) 
28 
Perhaps it will be useful to investigate here a 
little more completely the modifications that 
such a statement of problem imposes on the 
critical state. 

The variations of öz, of the lower interface 
always remains arbitrary so that the velocity 
u, must be zero because of the expression of 
the variation of the thrust (section 15). The 
general expression (10) of the velocity shows 
that the position of the maximum of the 
velocity is given by z = —«. Taking into 
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account that u = o for z = 2, and integrating 
this equation we obtain 


and, as it was observed already, z, is extremely 
small and the free surface can be assumed to 
[0] 
coincide with z = 0 so that / udz = qo. 
<0 
This condition can be immediately written 
in dimensionless form. Calling h the height 


of the current (h = — z,) and // the parameter 
(Ao/e) (gh°/q°), 

I I & 

LA > 6 Bi 


We must observe that up to now we have 
not used the surface condition; if we designate 
by m = «/h the relative depth of the velocity 
maximum, we will have in the most general 


manner and for all surface conditions 


Taine Stains (14) 


In particular, if the maximum is at the surface 
we find for the critical parameter the value 
II = 6 previously computed. 

Let us consider now Eq. (11). Noting that 
the velocity is zero for z,, we have 


O = 008% — Pin 
and this relation becomes 
QU? = p + 092 — ga (00 — 0) 
In particular, on the free surface where pro 
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Equation 12 becomes, by introducing the 
preceding value of 2, 
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and in its dimensionless form 


48 I (; Ao 
= m 
2 8) Me ONG 0 
But the velocity at the surface is given by 


Eq. (13) in which we can again assume that 
2, = 0; thus, 


u = À 
4h 


(1—2 m) 

We can note here (an obvious property of the 
parabola) that if m > 1/2 the velocity at the 
surface becomes negative; we have to exclude 
such solutions from a physical viewpoint, 
because they correspond to an instable gradient 
of the density and the appearance of densities 
foreign to the problem; m can then vary 
from — oo to 1/2, but beyond that, the mathe- 
matical solution has no physical significance. 

Writing now the imposed condition 


u2 I (; Ao 
= 10 | = 
Sehe 3 Xk D RD 


and by calling X the ratio od/hAo we get 


TT? dor 


2 
ae (1 — 2 m) 


The elimination of m between Eqs. (14) and 
(15) leads to 


We must note, however, that the auxiliary 
quantity X contains the unknown h sought. 
We shall separate the given quantity and the 
unknown by setting 


which gives 


We have then in final form 


ana td Man ee I 
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to which we must add the relation 


as giving a clear picture of the velocity distri- 
butions in relation to the relative depth of the 
maximum. 

As K is not very expressive, we can write 
it also in terms of the height and velocity of 
an initial state of uniform velocity and con- 
stant density gradient. It can easily be shown 
that in this case z, = (1/2) (Ao/o) ho and as 
u,?/2 ¢ = q0?/2g hu’, we obtain by setting 
II, = (A@/e) (g ho®/qo?) the form 


Le I 
sen aa) 


Thus, the critical value of IT depends this 
time on //,, and one could interpret the pre- 
ceding results as a tendency of an initial state 
II, toward a (corresponding) critical state 
characterized by 11. 

The study of the case where the level of 
the free surface is imposed 


Ro 


zı=d 


leads to similar results. If we introduce again 
a parameter K = V/(o/Ao)? (gd?/go?) the cri- 


tical value of JT is given by 


K= VA (= +2) 


1140003 


and K is related to a critical state of constant 
velocity and constant density gradient by 


i 8/77 
> 5 VI 


Summary and perspective 


As has been shown in the foregoing pages, 
the momentum principle and the energy 
principle of mechanics permit two different 
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integrals to be written for the cross section of 
a gravity current of constant width. Because 
of external shear on the one hand and energy 
loss on the other, both integrals must decrease 
in magnitude in the direction of flow. More- 
over, if we consider all the possible states 


which the steady flow of such a current can . 


take, each of the integrals is found to have a 
minimum magnitude corresponding to a 
certain critical thickness of the current and a 
certain critical velocity distribution. The 
critical states indicated by the two principles 
are not identical, although qualitatively com- 
parable. 

The expression “all the possible states” 
requires particular attention in the case of 
“ surface currents. If some condition restricts 
the vertical displacement at the free surface, 
the critical state will be modified by this 
condition. In general it will result in a critical 
velocity distribution in which the maximum 
is reached at some distance from the surface. 
If this condition does not exist, however, the 
maximum will occur at the surface. 

The reader might well ask what precise 
physical significance these critical regimes have. 
They appear to play the same role as a road 
sign on a highway, in that they provide a 
sense of direction — a prediction as to a 
future stage — but unfortunately without any 
indication as to the distance which remains 
to be traversed before that stage is reached. 
In the early stage ‚however, it is perhaps, the 
local effects of the surroundings which actually 
govern the velocity distribution rather than 
the perception of a low-ebb limit beyond 
which further travel is impossible. In other 
words, the critical stage may be only a hypo- 
thetical terminus — at least for the assumed 
nature of the current (very small curvature, 
hydrostatic pressure distribution, etc.). 


In open-channel flow, nevertheless, the 
critical stage is generally reached at such a 
catastrophic change in boundary conditions 
as an abrupt drop, and the effect of this change 
continues upstream in the form of a back- 
water curve. Not only may there well be 
analogous catastrophes in the case of gravity 
currents, but the critical regime may then 
also be an element of more general “‘back- 
water” effects. This hypothesis remains to be 
elaborated, just as we still have to verify the 
existing theoretical results through additional 
experimental evidence from nature and from 
the laboratory. 

A further question as to marine currents 
arises in connection with their freedom of 
lateral change as well as vertical. The answer 
to how and why lateral equilibrium is estab- 
lished must await further progress in the 
theory. Perhaps one will be tempted to con- 
clude from these remarks that what is known 
is little in comparison with what remains to 
be learned — but is this not the general situa- 
tion in all science? 
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Determination of the Lateral Eddy Diffusivity in the 
Climatological Mean Gulf Stream 


By H. STOMMEL, Woods Hole Oceanographic Institution 


(Manuscript received 23 December 1950) 


The spread of positions of the Gulf Stream as 
exhibited in fig. 1 of the paper by Fuctister 
and WORTHINGTON (1951) published in this 
issue of Tellus suggests that one should be able 
to compute a rough value of the lateral eddy 
diffusivity produced by the meandering and 
shifts of position of the stream, although the 
number of points at any section, taken across 
the stream at any distance along the axis of the 
stream is quite small for statistical purposes. 

If o represents the standard deviation of the 
position of the actual stream from the mean 
axis of the stream at any point along the axis, 
then, roughly speaking, the eddy A should be 
given by the equation 


A=o V/2L, 


where L is the distance of the section from 
Cape Hatteras and V is the mean velocity of 
the stream. o has been evaluated for two 
section for which L is 500 kilometers and 1,000 


kilometers. The results of the computation 
are given in the following table. 
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It must be remembered that this determina- 
tion of the lateral eddy diffusivity is based on 
the assumption that the actual Gulf Stream 
is a turbulent flow. The numerical values are 
roughly of the same order of magnitude as 
those used by Munk (1950), though somewhat 
smaller. 

The result of this computation seems to 
support the contention that the western in- 
tensification exhibited in the Munk-Stommel 
theories are climatological or statistical in 
nature and that the very thin, filament-like 
nature of the individual turbulent element of 
the Gulf Stream may be due to another cause. 


REE EE RIE INS ES 


FUGLISTER and WORTHINGTON 1951: Some Results of 
a Multiple Ship Survey of the Gulf Stream. Tel- 


IUT ETC SD Ee 


Munk, W. H. 1950: On the Wind-driven Ocean 
Circulation. Journ. of Met, 7. pp. 79—93. 


1948—1949 Trials of the Schaefer-Langmuir Cloud-Seeding 
Technique in Hawaï” 


LUNA B. LEOPOLD and WENDELL A. MORDY 


Abstract 


Fifteen tests during 1948—49 involving thirty-seven individual cloud inoculations over the 
islands of Lanai and Molokai are described. With regard to the rôle of the dry ice in causing the 
observed rainfall, the tests are too few in number to be considered conclusive, but the large rain- 
fall values associated with seeded clouds are of particular interest. The largest rains seem to be 
associated with clouds where cloud top temperatures are slightly colder than o° C and such 
clouds are also relatively thick. The relative importance of cloud thickness and temperature 


cannot be determined from the available data. 


There exists at the time of this writing 
divergence of opinion concerning the efficacy 
of seeding clouds with dry ice for the produc- 
tion of rain. 

The current status of the problem was 
succinctly stated by BERGERON (1949). “The 
question... is... no longer whether pre- 
cipitation can be released artificially, but rather 
whether such a release can produce an appreci- 
able amount of rainfall, and when and where 
this could be done.” 

The early experiments in Hawaii have 
already been reported (LEopotp and HALSTEAD 
1948) and the present paper presents the results 
of further seeding experiments conducted in 
a more or less comparable manner. The 
consistent location of seeded clouds from day 
to day, the large number of rain gages, and 
high elevation of the freezing level make the 
Hawaiian experiments of some interest, though 
both the earlier and the current sets of data 
suffer the disadvantage of no observations 
with radar. 

In all but one of the trials discussed here, the 
. dry ice was ground to a size resembling granu- 


ı Published with the approval of the Director, as 
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lated sugar. The plane was a twin-engined 
Beechcraft in which 400—s00 pounds of dry 
ice were carried on a given flight. The dry ice 
was dispersed from a hopper having a sliding 
door, and from which 150 pounds of the 
granulated material could be dropped out in 
a time period of five to seven seconds. Flying 
at a speed of 140 mph the seeding rate was 
approximately 12.5 pounds per 100 feet of 
horizontal flight. The cloud towers seeded 
were ordinarily of such a width horizontally, 
about one-quarter mile, that the six seconds 
during which the dry ice was dropped was 
approximately the flight time across or through 
the top of the tower. For a seeding run the 
plane was flown just through or just over the 
cloud top. 

The clouds seeded were those which occur 
over the central part of the Island of Lanai or 
over the Maunaloa area of the Western part 
of the Island of Molokai. These clouds on test 
days were cumulus congestus associated with 
the cloud line formed by the interaction of 
tradewind and sea breeze. These clouds have 
been described in detail in another paper 
(LeoroLD 1949). The cloud type is important 
in considering the results because the con- 
vergence in the lower levels where the op- 
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posing winds meet provides a source of ver- 
tical motion in the upper cloud which is not 
found in the common air mass cumuli. It is 
possible, then, that these clouds possess the 
prolonged “path of fall” of droplets referred 
to by BERGERON (1949, p. 34), even though 
the total cloud thickness seldom exceeded 
12,000 feet at the time of seeding. 

It can be stated with some certainty that 
when the cumulus clouds are capped by a tem- 
perature inversion of 1° C or more as is the 
normal case in Hawaii, seeding with dry ice 
either produces a partial dissipation of the 
cloud or no noticeable effect; the clouds do 
not increase in height and the rain, if any, is 
of insignificant amount. This is based on a 
large number of observations during the 
1947—1948 tests including not only tests made 
by the authors but many others run independ- 
ently by various companies in Hawaii. The 
conclusion was sufficiently well established 
_during the first year that all later tests reported 
here were conducted on days when the Hono- 
lulu 1500 Z (osoo LST) radiosonde indicated 
the absence of the tradewind subsidence in- 
version or that the inversion was less than 1° C. 
Changes in the strength of the subsidence in- 
version occur rapidly in the Hawaiian area. 
On several occasions the osoo LST raob 
showed no inversion, on the basis of which a 
seeding flight was made. However, by 1000LST, 
after the flight had already begun a definite 
inversion had formed, as could be seen by the 
concordance of cloud tops and by temperature 
measurements taken in flight. 

Our experience, furthermore, indicates that 
large horizontal wind shear through the levels 
between base and top of cloud tends to blow off 
the top of the cloud, and rain will not result. 

An initial depression of the cloud top im- 
mediately following seeding has been observed 
by nearly all experimenters. Experience in 
Hawaii indicates that the cloud will recover 
from this initial depression and will grow 
even higher than its pre-seeding condition 
only if the cloud was actively growing in 
height at the time of seeding. Growing clouds 
over Hawaii usually have “harder”, more 


well-defined tops than those not in the process, 


of growing which tend to be wispy-edged. 
The wispy edges, the authors believe to be 
either a dissipating phenomenon or a result of 
ice-crystal formation in the clouds. 


It is typical for moisture values to be moder- 
ately large below the subsidence inversion, 
decreasing sharply above. On days when the 
inversion is absent, moisture occurs through a 
much thicker layer, and experience points 
toward large values through a deep layer as a 
feature conducive to rain from a seeded cloud. 

The Hawaiian experiments share with many 
others the difficulty of defining a practical and 
definite control against which the results of 
seeding can be compared. The tests reported 
here were all made on days of small or no 
temperature inversion. On such days, all 
clouds in the area tend to be thicker than usual. 
The same conditions promote the growth of 
the sea-breeze clouds. For these reasons, it is 
obvious that the days chosen for test are the 
days on which natural rainfall is most likely to 
occur on the test islands, Molokai and Lanai. 

On the other hand, it is a fact widely known 
in Hawaii that even on such days, when 
unusual build-up of cumulus clouds can be 
seen everywhere, dark threatening sea-breeze 
clouds remain nearly all day over the Islands of 
Molokai, Lanai and certain other areas, but 
seldom does more than a light shower occur. 
One reason for this undoubtedly lies in the 
fact that such a cloud, whose appearance 
suggests imminent rain, often extends nearly 
to the freezing level but seldom into it. The 
freezing level over Hawaii usually lies between 
14,000 and 18,000 feet, while the cloud tops 
on days of no inversion generally lie between 
10,000 and 13,000 feet. On such a day a few 
cumulonimbus usually can be seen somewhere 
in the Territory. It appears, then, that these 
conditions are unusually favorable for the use 
of dry ice as an agent to start the rain process 
which is potentially near but whose unaided 
initiation is relatively infrequent. 


Result of tests 


In the 1948—1949 series, tests made on 6 
days over the Island of Lanai and 9 days over 
the Island of Molokai include adequate data. 
Certain additional tests are omitted because 
necessary data are not available or because 
circumstances made the conditions of test 
sufficiently different from the standard that 
the tests cannot be considered comparable. 


For example, some clouds were seeded which 


lay entirely over the ocean and though rain 
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was observed, it could not be measured. In 
certain others, wind and temperature condi- 
tions changed so rapidly during the test that 
the Honolulu upper-air data were considered 
completely unrepresentative of the test area. 

On a given day of test, one to five seeding 
runs were made. Data on individual runs, in- 
cluding summaries of cloud changes observed, 
are included in Table 3 appended to this report. 
The rainfall quantities measured are sum- 
marized in Table 1. These data represent the 
arithmetic mean rainfall in five gages except 
where indicated. 


Table 1. Summary of all Hawaiian experiments, 
1947-1949. 
Number of Test Days: 27. 


Cases of Cloud Tops Cases of Cloud Tops 


Warmer than o °C Colder than o °C 


No. Rain at Ground No. 
of or of 


Rain at Ground 
or 


@ases Other Notes Cases Other Notes 
I 1.251 I | 2.182 
1 | 0.03 | I 1.952 
OMIRTEICe 1.88 M 
0.78 L 
2 | Virga I 0.602 
3 | Cloud dissipated I 0.50 
6 | No change observed I 0.25? 
19 ip Trace 
1 | Cloud dissipated 
8 


1 Cloud went above freezing level after seeding and 
during rain. 

2 Mean of 5 gages in seeded area; rainfall figures 
without superscript are for an individual gage under 
seeded cloud. 

M Molokai; L Lanai. 


It will be noted that the rainfall quantities 
measured in association with the experiments 
reported here are considerably larger than 
those reported by Orr et al. (1950), Coons 
et al. (1948 a) and most other workers. In the 


Lanai test of June 28, 1948, 1.75 inches fell in’ 


45 minutes at Lanai City at the outer edge of 
the cloud. A total of 3.40 inches fell in the 
Palawai Basin directly beneath the seeded 
cloud. This rain began 19 minutes following 
seeding. It is the largest rainfall on any day in 
June in a 45-year record on that island. ; 


1 Gages used in this analysis: Lanai gage numbers 
671, 672, 674, 681, 683; Molokai gage numbers 513, 
$14, $18, $19, 522. Locations indicated by maps of Lro- 
POLD, BURN, STIDD (1948). 
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The interpretation of the data shown in the 
table will be facilitated by a discussion of 
certain problems which will immediately occur 
to the reader. 


Did any rain precede seeding, and was 
rain occurring in the general area? 


This question was: stressed by Coons, 
Gentry and GUNN (1948 c) who noted in a 
particular test that radar echoes from snow 
were existing or were forming in the imme- 
diate vicinity at the time of seeding. ” Without 
the facilities utilized ..., especially the radar, 
it probably would have been concluded that 
the one inch of snow was definitely a result of 
the cloud-seeding operation (p. 18).” 

Theirs is a cogent argument. From the 
standpoint of agriculture, however, it is not 
necessary to prove that dry ice can precipitate 
rain when no rain would have fallen. The 
question is rather: was the amount or timing 
of rainfall affected by seeding? 

The Hawaiian experiments were all run on 
days of unusually high probability of natural 
rainfall, and rain usually did occur in the area 
near or at the time of rain from the seeded 
cloud. However, rainfall amounts under the 
test clouds were generally larger than those 
in the vicinity. 

As mentioned earlier, it is difficult to 
define a practical and definite control by which 
we may compare the results of seeding against 
circumstances which might have occurred had 
no seeding taken place. As a step in this 
direction, results for the 8 experiments where 
clouds colder than freezing were seeded are 
presented in Table 2 and compared with rain- 


Table 2. Mean rainfall amounts measured by 
specific rain gages on days of test and similar 
days of no test. 


00% 


OI—.09 

.10—.99 
> 1.00 
Total 


Rainfall 


(Rainfall values are average 
of $ gages, Lanai) 


Non-seeded days: 


NOMO days ne artecras 26 25 Die Sel 77 

ey ore WOVEN TEN Re 34 32 27 7 
Seeded Days: 

INO ROMANS RER 3 2 I 2 8 


0/ 


of total days......] 
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Fig. 1. Rainfall amounts associated with seeded clouds of given temperature and thickness. Temperature given 


are at the level of CO: release, approximately at the top of the cloud. 


fall amounts occurring on days where no 
seeding was done and similar upper-air and 
cloud conditions existed. 

Because there are only 8 such cases, a sta- 
tistical test for significant differences in the 
rainfall frequency distributions on these days 
versus non-seeded days will not yield a sig- 
nificant result. The results are nevertheless of 
interest even if definite conclusions cannot be 


reached. 


In what ways did rain occurring after 
seeding appear to be related to seeding? 


As previously mentioned, rainfall in each 
case followed the seeding by only a short 
interval of time or rain intensity increased. In 
most instances, rain endured for less than an 
hour. In the case of Test 1 L and 5 L, however, 
rainfall continued for approximately three 
hours. The heaviest rainfall intensity in both 
of these instances occurred between 15 and 30 
minutes following seeding. Clouds having 
hard or well-defined edges at the top were 
chosen for seeding with the idea that rain in 
such clouds had not yet developed. Lacking 
radar this was the best evidence concerning 
the state of the moisture droplets in the cloud. 
Other experimenters equipped with radar 


equipment have indicated that this feature is 
indicative of the lack of formation of ice 
crystals. within the cloud (Orr et al. 1950). 


What Was the Relation between Rain Cloud 
and o° C Isotherm? 


Data plotted in fig. 1 indicate that for the 
15 days of experiments in 1948—1949 re- 
ported here, a difference in the rainfall amounts 
exists between those where the cloud-top 
temperature was warmer than freezing and 
those where the cloud-top temperature was 
colder than freezing. Since the freezing level 
in Hawaii tends to be at roughly a uniform 
height, this represents, as indicated in fig. 1, 
not only colder temperatures but also great- 
er cloud thickness. A correlation of cloud 
thickness with observed rainfall was indi- 
cated by the previous Hawaiian experiments 
(Leoporp and HALSTEAD 1948), and its 
importance has been stressed by VIERHOUT 
(1950). Since even the coldest temperatures 
observed at the cloud tops were only — 5° C, 
it is difficult to evaluate the relative importance 
of temperature and cloud thickness. 


Conclusions 


No definite conclusions can be drawn from 
these data owing to the small number of 
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Table 3a. CO, Tests 


Island 
Data from 1500 Z Sounding Honolulu Height and Thickness 
Inversion 
Nee as Ps asia Pena AUS ETS NE a es of Quantity 
5 : rop LST of CO, 
Pressure Strength! & Altitude Bas Top 
At Base at Moisture 
Break 
(Began 3 min. 
scatter 
1 M 6/28/48 701 mb 0.4 700 mb 135531 ESE) 100 Ibs 3,000 ft 15,500 ft 
x 9,800 ft 9,850 ft 14:03 100 3,000 15,500 
2M 6/29/48 730 mb 0.3 730 mb I2: 00 150 3,000 8,000 
8,775 ft 8,775 ft I4: 20 Iso , 3,000 8,000 
3 M 7/9/48 680 mb 14) 700 mb 13 07 100 2,900 11,000 
10,600 ft 10,500 ft 13: 30 100 2,900 12,100 
13:36 so 2,900 11,300 
15: 34 100 2,900 8,000 
15:39 so 2,900 8,300 
4M 8/3/48 632 mb 0.8 755 mb 14: 22 100 2,500 9,500 
12,450 ft 7,900 ft 14:32 100 2,500 10,000 
14:45 100 2,500 10,500 
5 M 9/22/48 659 mb 1.4 631 mb 13402 100 4,000 11,500 
11,425 ft 12,950 ft FICO 100 4,000 12,500 
| 13:29 100 4,000 13,000 
6M 9/24/48 525 mb 0.7 525 mb. 12,23% 75 2,500 13,800 
17,100 ft 17,100 ft 12:47 75 2,500 14,500 
TRS Iso 2,500 14,800 
7M 11/9/48 608 mb 0.7 610 mb 134.19 150 2,000 8,160 
13,425 ft 13,375 ft 
8M 11/27/48 850 mb 0.3 None 14: 33° 150 2,000 14,000 
None 15: 048 Iso 2,000 13,100 
9 M 6/13/49 585 mb 0.5 573 mb PRESS 100 2,500 12,900 
14,400 ft 14,900 ft 


0, Top — O; Base _ 
4 


1 Figures given are approximately 


2 Fine granules. 
3 Coarse granules. 
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in 1948—1949. 


of Molokai. 
Rot tage Alt. & Temp. at 
Level Ice Release 
wert Mean 
Temp. mA  OnunGrlonuds sa midge Rea wn fa 101 Rainfall 
Thickness} Alt from 5 Gagess 
; Honolulu 
Sonding 

12,500 ft} 14,300 ft] —2.0°C| 12: 00 LST weather report from Molokai indicated rain showers. Rain] 1.59 in. 

12,500 14,000 — 1.2 was not visible to pilot at altitude. Ice formed on plane while seeding. 
Heavy rain was seen from plane at 14: 33 LST under seeded towers. 

$,000 8,000 + 0.0 No rain. 0.00 

$,000 8,000 == 0:0 

8,100 10,700 + $.5 Light showers began at 13: 10 LST from cloud seeded at 13: 17 LST.| 0.23 

9,200 12,100 + 4.0 A sudden increase in rain could be seen at 13: 18 LST by observers on 

8,400 II,500 + 5.2 Lanai (shown in photographs). At 13: 47 LST a further increase of 

5,100 8,000 + 8.0 rain was noted. At 13: 56 LST rain covered most of West Molokai 

5,400 8,300 + 8.2 as seen from Lanai. 

7,000 9,500 + 5.5 Cloud dissipated in 2 minutes after final seeding. No rain associated] 0.00 

7,500 10,000 + 4.5 with seeded cloud. Clouds in vicinity showed wind shear with tops 

8,000 10,500 + 3:0 and bases separating. 

7,200 II,300 + 2.1 No rain. Cloud dissipated. 0.00 

8,200 12,100 + 2.0 

8,700 12,400 + 2:0 
11,300 13,500 — 1.1 Rain was falling from cloud at time of seeding. At time coincident with] 0.25 
12,000 14,000 — 2.0 first seeding, observer on Lanai noted darkening of cloud and a second 
12,300 14,300 — 2.5 shaft of rain from base of cloud. At 12: 43 LST rain covered 10 times 

the area prior to seeding. Rain activity decreased at 12: $5 LST and a \ 
new shower began at 13: 04 LST. Cloud tops flattened out beginning 
ap 132 O4 EST. 
6,160 8,000 + 8.0 Cloud dissipated rapidly. 0.00 
12,000 14,500 — 3.8 Rain at time of seeding. Cloud tops lowered after first seeding (14: 40] 0.00 
11,100 12,800 == où LST); at rs: 05 LST rapid dissipation of cloud could be seen from 
Lanai. 

10,400 12,500 + 0.2 Showers over seeded area at 12: 15 LST. Of short duration and small} 0.60 
amount. Heavy rain observed from plane under seeded cloud at 
1 yi) SIE 


4— 006782 


so 


No. & Date 


1 L 6/28/48 


2L 11/9/48 


3 L 11/18/48 


4 L 11/26/48 


5L 1/10/49 


6L 6/13/49 
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Table 3b. CO, Tests 


Island 


Data from 1500 Z Sounding Honolulu 


Pressure 
at Base 


701 mb 
9,800 ft 


608 mb 
13,425 ft 


550 mb 
15,900 ft 


691 mb 
10,200 ft 


635 mb 
12,350 ft 


585 mb 
14,400 ft 


Inversion 
Pressure 
& Altitude 
1 
Strength at Moisture 
Break 
0.4 700 mb 
9,850 ft 
0.7 610 mb 
13,375 ft 
0.2 soo mb 
18,300 ft 
0.3 692 mb 
10,150 ft 
0.5 735 mb 
8,600 ft 
0.5 573 mb 
14,900 ft 


9; thee — 0; Base R 


= legen, Falikaal cha sap pro xninate |) 


? Ice was scattered more slowly than for other tests. 


Times of 
Drop LST 


12:56 LST 
73073 
13: 19 


12255 
13: 01 


Quantity 
of CO, 


100 lbs 
100 
50 


100 


100 
100 


Height and Thickness 


Base 


3,000 ft 
3,000 
3,000 


1,800 


2,000 
2,000 
2,000 


unknown 
unknown 


3,000 
3,000 
3,000 
3,000 


3,000 
3,000 


Top 


15,500 ft 
14,500 
14,500 


12,000 


10,500 
10,600 
10,500 


7,000 
6,800 


14,000 
14,000 
12,000 
12,000 


14,400 
14,400 
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in 1948—1949. 


of Lanai. 
Alt. & Temp. at 
f : 
of Cloud Level Ice Release 
— Mean 
Temp. Reeması son Cloud ATTAMRR En ERaNNl Rainfall 
Thickness Alt. from 5 Gages 
Honolulu 
Sounding 
12,500 ft} 13,800ft | — 1.0°C| Between 1 and 2 in. fell previous day. Cloud built rapidly in 2! 1.95 in. 
11,500 14,000 — 1.2 min. following ice but 14 min. after seeding, top of cloud was 
11,500 14,000 1,000 ft. lower. Kain at Lanai City began at 13:15 LST and 
in Palawai Basin (seeded area) about 13: ı5 LST. 
10,200 12,000 + 1.2 Rain had occurred in hour previous to test. Tops of CU were) 0.06 


flattened into As at 7,000 ft. At 13: 53 LST plane flew under 
seeded arca and encountered light rain and moderate turbulence. 


Cloud top dropped 500 ft 4 minutes after 14:27 LST. Large} 0.03 


8,500 10,500 35 
8,600 10,600 + 3.2 central portion of cloud dissipated. At 15:00 LST rain was ob- 
8,500 10,000 + 3-5 served from plane and ground directly below seeded cloud. No 
other rain on Lanai or Molokai. Stopped raining at 15: 30 LST. 
6,900 + 9.8 No rain. 0.00 
6,700 10.2 
11,000 14,000 — 5.0 Small showers fell in morning but were not associated with CO,| 2.18 
11,000 14,000 |— 5.0 drop in A. M. Cloud seeded was not highest cloud in area but 
9,000 12,000 |— 3.0 was in formative process with well defined hard tops. Cloud 
9,000 12,000 3.0 was small in area. Rain began as a single shaft directly below 
point where plane was seen to emerge from cloud (700 ft below 
top of cloud). Rain rapidly spread over wide area. 
11,400 14,400 |— 3.2 No indication of rising cloud tops after seeding. Top of cloud] o.12 
11,400 14,400 | — 3.2 dissipated after second seeding and became separated from lower 
portion. Tops became 8,000 ft by 13:25 LST. Moderate rain 
was observed at 13:30 LST when plane flew under seeded cloud. 


42 LUNA B. LEOPOLD and WENDELL A. MORDY . 


tests and the lack of suitable control. Never- 
theless, it is important that all pertinent data on 
the subject of cloud inoculation be added to 
the literature. The present experiments are of 
particular interest because of the large rainfalls 
which were experienced under the seeded 
clouds and very possibly related to the in- 
oculation. 
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Note on the Dynamics of the Antarctic Circumpolar Current’ 


By W. H. MUNK? and E. PALMEN? 


Abstract 


Unlike all other major ocean currents, the Antarctic Circumpolar Current probably 
does not have sufficient frictional stress applied at its lateral boundaries to balance the 
wind stress. The balancing stress is probably applied at the bottom, largely where the 
major submarine ridges lie in the path of the current. The meridional circulation pro- 
vides a mechanism for extending the current to a large enough depth to make this 


possible. 


Under the assumption that the wind stress on 
the sea surface is balanced by the frictional stress 
against the sides of the ocean basins one can arrive 
at a fairly satisfactory picture of the large scale 
circulation in various oceans (MUNK, 1950; MUNK 
and CARRIER, 1950). The computed transports differ 
at most by a factor of two from the observed 
transports. In the case of the Antarctic Circumpolar 
(AC) current a similar reasoning leads however to 
a computed transport one hundred times the ob- 
served transport. 

In its simplest form we may regard the AC 
current as an eastward flow on a plane tangent to 
the earth at the South Pole. The flow is induced 
by constant eastward winds, and depends only on 
the distance r from the pole on this plane. In this 
system a balance between the wind stress t and the 
lateral friction is expressed by 


d (+ =) 
t+A - - =|) 26 


dr \ dr r 


where A is the lateral kinematic eddy viscosity and 
M the eastward mass transport across a normal 


1 Contribution from the Scripps Institution of 
Oceanography, New Series No. 515. This work has 
been supported by a contract with the Office of 
Naval Research. 

2 Institute of Geophysics and Scripps Institution of 
Oceanography, University of California. 

3 Academy of Finland. 


vertical plane of unit width extending from surface 
to bottom.‘ The solution which vanishes at the 
Antarctic continent (r =o) and at some other 
latitude (r = rı) is 


T £ x Gr hoe r, 
= = 172 To? : = In . 
18 A (ee Syn If ity 


Setting 7 = 2 dynesiem 2, A = 108 cm sec=4 and 
placing the boundaries at 70° S and 45° S latitudes 
gives l= 57% 100 otsee= At the narrowest sec- 
tion, Drake Passage, the boundaries are at 65° S 
and $5 7's latitudes, and 1 = 1625" 8 seche 
observed transport is 1014 g sec~'. The discrep- 
ancy is not materially altered by more elaborate 
calculations involving spherical coordinates and 
allowing for a variation of wind with latitude. 


4 Strictly speaking we should replace the relative 
transport M by the absolute transport M + wrh, where 
© is the earth’s angular velocity and h the ocean 
depth. However, the second term which represents 
solid rotation drops out upon differentiation. 
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Among the possible explanations for the discrep- 
ancy are: 
1. The value of A pertaining to the AC may have 
been taken too small, and should possibly be one 
to two orders of magnitude larger than 108 cm?sec~!. 
Various methods for other large scale ocean currents 
have led to values from 107 to 10° cm? sec-r. 
2. The wind stress is balanced by stress against the 
bottom, rather than against the lateral boundaries 
as has been assumed. Somewhat artificially we may 
divide this bottom effect into (2a) skin friction, and 
(2b) the mountain effect. 

2a. A bottom velocity so great that the wind 
stress is balanced entirely by skin friction, is given 
by Sverprup et al. (1942, p. 480), 


T Z + Zo 
EN, \/ = In ( ) 
0 zo 


Setting t = 2 dynes cm, 9 = 1.0 g cm, 2 = 2 
em, and == 100 cm gives 0, = 14 cm sec at 
100 cm above bottom. This is approximately the 
maximum surface velocity. The value for the 
“roughness length” zo is based on unpublished 
measurements by Revelle and Fleming in San Diego 
Harbor. But this process could at best lead to a line- 
ar velocity gradient between top and bottom of 
Dico sec 4 km ==03.5 105 seez" In order 
for the resulting vertical stress, 3.5 + 1075 A,, to 
equal a wind stress of 2 dynes cm”? at the sur- 
face the vertical eddy viscosity A, would have 
to equal 6- 104 g cmt sec~!. Reported values 
are generally less than 2 - 103 g cm! sect in 
the upper few hundred meters, and even less at 
greater depth (SVERDRUP et al. 1942, p. 482). 
2b. The foregoing equation refers to the statis- 
tical effect associated with “small” irregularities in 
the sea bottom. In meteorology, the mountain 
effect associated with the few largest irregularities 
on the earth’s surface has been found at least as 
effective as skin friction (WIDGER, 1949), but in 
oceanography the mountain effect has been neg- 
lected. Let dh/dx denote the eastward (down current) 
slope of the sea bottom. The eastward component 
of force exerted on an element of bottom of area 


ds equals 


py (dh/dx) ds, 


where p, is the pressure along the bottom. For 
small slopes, ds can be taken as the area projected 
on a horizontal plane. The net effect around a 
circle of latitude of circumference C vanishes if the 
pressure is constant (/ dh = 0), or if the depth is 


q 


constant (dh/dx = 0). If however the bottom pres- 
sure on the up current side exceeds on the average 
the pressure at the same level on the lee side by 
Ap, then the stress has a mean value of the order of 


I 
— Ap, Al 
C Pp An 


which could balance the mean wind stress. Here 
Ah is the sum of the heights of the sea mounts. At 
60° S the current encounters the Kerguelen ridge 
(3 km), the Macquarie ridge (2 km), the South 
Pacific ridge (1 km), and the South Antillean arc 
(4 km). The currents must diminish over each 
ridge and accelerate between ridges, but on the 
average the mountain effect may balance the wind 
stress. Setting 4h = 10 km, C = 18,000 km, 
Tt = 2 dynes cm? gives À p, © 4,000 dynes cm? 
= 4 dynamic centimeters. At the surface the hori- 
zontal variations in pressure are of the order of a 
dynamic meter, with the larger values occurring on 
the up current side. The wind stress can therefore 
be accounted for by horizontal pressure gradients 
at the bottom similar in phase to those observed 
near the surface, and amounting to a few percent 
of the surface gradients. 

If the wind stress is to be balanced by stresses 
along the sea bottom, the AC current must be deep 
enough to reach to the bottom. Sverdrup has pointed 
out that the sharp leftward displacement of the 
surface current over the submarine ridges would 
indicate that this is the case (SVERDRUP et al, 1942, 
pp. 466—469). We may consider two possible 
explanations: 

1. Each layer induces, by turbulent interchange 
of momentum, motion in the layer beneath, and 
in this manner the wind stress is transmitted to the 
sea bottom. Because the fetch of the wind is essen- 
tially infinite this mechanism is more effective in 
the AC current than elsewhere. 

2. The meridional circulation exports absolute 
angular momentum in the upper layer, and imports 
it in a lower layer. Suppose the loss of angular mo- 
mentum in the upper layer alone associated with 
a flow of Q g sec”! northward across the 45° S 
parallel of latitude is balanced entirely by the wind 
torque over the ocean between the Antarctic con- 
tinent (70° S) and 45° S: 


Q (© R cos 45°) (R cos 45°) = 
FOÈRS 


= ft (27 R? cos y) (R cos y) d y, 
45°S 


where & is the earth’s angular velocity and R its 
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radius. Setting t= 2 dynes cm? leads to Q = 
3 X 103g sec". This value is in agreement with the 
observed northward transport in the upper $00— 
1,000 meters, due mostly to the Benguela and Peru 
Currents. If we suppose for the moment that the 
difference between precipitation and evaporation 
over the whole region south of the latitude 45° S 
and in addition a possible decrease of the Antarctic 
ice cap could provide all this water, then the wind 
stress could be accounted for by the export of ab- 
solute angular momentum into the southern oceans. 
Actually not more than about 0.3—0.5 per cent of the 
above value of Q could be provided in this ways, 
so that virtually all of the water lost in the upper 
layers must be returned to the Antarctic Sea at a 
lower level‘, apparently at depths between 1,000 
and 4,000 meters. Conservation of angular momen- 
tum then requires an increasing cyclonic (eastward) 
circulation at these great depths with decreasing 
distance from the pole. Deacon’s analysis of the 
DISCOVERY observations in fact shows a re- 
markably uniform horizontal density gradient 
throughout most of the water column, whereas, 
for example, Iselin’s sections across the Gulf Stream 
reveal a sharp decrease in horizontal gradients at a 
thousand meters (compare figs. 159 and 186, 
SVERDRUP et al. 1942). From this point of view the 
principal action of the wind is not so much to induce 
cyclonic circulation throughout the water column, 
as to prevent anticyclonic circulation of the north- 
ward drifting surface waters. 

To summarize: deep water in the Antarctic ocean 
drifting towards the axis of the rotating earth 
develops an easterly (relative) flow. The retarding 
pressure of the submarine ridges against the deep 
current balances the stress exerted by the wind on 
the water’s surface. 


5 Provided a reported secular rise in world-wide 
sea level by I mm per year is due alone to melting 
in the Antarctic. 

6 An additional complication is the northward flow 
along the bottom of the heavy water formed near 
the Antarctic continent during freezing, and account- 
ing for most of the bottom water in deep basins the 


world over. 


The foregoing speculation regarding the dyna- 
mics of the AC current would make it differ radi- 
cally from other major ocean current systems, which 
are limited essentially to the upper 1,000 meters, and 
for which the stresses against the lateral boundaries 
are adequate to balance the wind stress. The differ- 
ence is due largely to boundaries which prevent the 
zonal flow of water around the earth. Assume that 
the Drake Passage should close, and that the Mac- 
quarie Ridge should emerge from the surface. Then 
the water in the confined sector would pile up 
against the downwind boundary. If the wind blows 
with equal strength at all latitudes, the resulting 
horizontal pressure gradient will balance the wind 
stress and there would be no circulation involving 
a horizontal gyre. If however the winds at higher 
latitudes are weaker or, as is the case, reverse in 
direction, the resulting wind torque would set up 
a gyre within this confined sector of the Antarctic 
Ocean. As an example, the Gulf Stream is the 
western component of a “subtropical gyre” main- 
tained by the easterly winds to the south and the 
westerly winds to the north. In a confined ocean 
the transport depends therefore on the curl of the 
wind stress (MUNK, 1950); in a zonal (Antarctic) 
ocean the transport has been shown here to be 
proportional to the wind stress itself. 

We have also shown the transport in a zonal 
ocean to vary as A-! provided the wind stress is 
balanced by lateral friction. In a confined ocean 
where such a balance does in fact exist the transport 
is independent of A (Munk, 1950)! This surprising 
result has to do with the west-east asymmetry in- 
duced by the earth’s rotation on a confined gyre. 
The center of the gyre (Sargasso Sea) is displaced 
to the west, and the western current (Gulf Stream) 
is pressed against the boundary. If A is very large, 
this asymmetry, and hence the shear against the 
western boundary, are small. If A is relatively small 
the asymmetry is very pronounced, and the shear 
very large. The product of A times the shear, i. e., 
the lateral friction, remains the same, and is so 
large in this western region that it almost balances 
the wind stress over the entire subtropical gyre. 
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Sporadic E Movements on 21 June 1949 


An interesting instance of sporadic E reflections 
occurred on 21 June 1949 over North America. 
These reflections, reported by radio amateurs 
operating in the 50 mc/s band during the hours of 
0000—0400 GMT, are portrayed in figures I—5. 
The various radio amateurs report whenever a radio 
contact is made by means of sporadic E. Such 
contacts are carefully screened to insure that propa- 
gation by means of the ground wave or the F2 
ionospheric layer is not involved. The midpoints 
of the great circle paths are then plotted as a dot 


in the case of single hoptransmission. With a double 
hop path, an “x” indicates the midpoint of each 
half of the great circle path. Each map represents 
the integrated number of contacts received dur- 
ing the one hour period centering on the hour 
indicated. 

The number of reporting amateurs during the 
entire period was about 300. It is seen that a rather 
extensive area of North America was blanketed by 
Es during the period shown. The greatest area 
covered was about 106 km?. No reports were re- 


Fig. 1. Hourly summary of sporadic E reported by oblique incidence during the hour centered on 0000 GMT. 


centered on oooı and 0002 GMT. 


Bic. 2 and 2. Hourly summary of sporadic E reported during the hour, 


5 


SE 


5 
Fig. 4 and s. Hourly summarv of sporadic E reported durine the hour centered on 0002 and 0001 GMT 
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Fig. 6. Summary map showing movement of spora- 
dic E during the period 0000—0004 GMT 21 June 
1949. 


ceived during the hours immediately prior to or 
after the time given. 

From the maps some type of drift motion at 
the altitude of Es (near 110 km) is evident. Whether 
such a drift is peculiar to this phenomenon (ie. 
solely an electron drift) or whether the reflection 
areas move about as part of a general atmospheric 
circulation pattern is not known with certainty. On 
the assumption that the general drift shown is an 
actual wind velocity, some indication of the average 
velocity of movement may be found. 

Two centers of Es are evident, both appearing at 
2400 GMT, 20 June 1949. The first group, B, 
centered over about southwestern Pennsylvania, 
moves from the ENE at an average speed of 240 
km/hr. Group À, first observed near central Mis- 
souri at 2400 GMT, drifts from the SE at an average 
speed of 250 km/hr. In both cases, the apparent 
velocity was greatest during the first hour (about 


350 km/hr) after which it decreased somewhat as 
dissipation and general disappearance of the “cloud” 
occurred. The average speed may be biased by the 
net number of radio observers operating at any 
particular time; however, in the above instance a 
sufficient number of amateurs were on the air both 
before and after the observational period to have 
reported Es had its effects been detectable at so mcys. 

It should also be noted that the term “average 
velocity”, which is determined from the movement 
of the center of mass of the plotted points, hardly 
describes the somewhat rapid thrust of the cloud A 
to the NW during the first 3—4 hours of observa- 
tion. The average velocity would be expected to 
be somewhat constant during the motion of the Es 
reflection area inasmuch as it is difficult to conceive 
of sudden movements or vortices in the atmos- 
pheric fluid at altitudes of 110 km. 

It is, of course, possible that an apparent motion 
only is displayed and that in reality the Es volume 
is dissipated at one side and produced at another, 
much like a banner cloud on the lee side of a moun- 
taintop. The mechanism for such an effect, however, 
is not immediately understood. A third hypothesis 
also exists: that the Es is formed by some sudden 
onslaught of particles from space, that these ionize 
the atmosphere in the vicinity of 110 km and that 
the ionized particles then move until dissipated. 
However, the maintenance of such an initially 
ionized volume for several hours under the known 
recombination rates of the E ionospheric regions is 
hard to explain. 

Appreciation is expressed to the radio amateurs 
in the United States for their wholehearted coopera- 
tion in supplying the observations, to Mr O. P. 
Ferrell in screening the data, and to Mr F. Crowley 
for plotting. 


N. C. Gerson 
AF Cambridge Research Laboratories. 
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